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ABSTRACT 


The equations of motion for a ducted fan hovering device 
are developed and programmed in a computer simulation. 
Experimental aerodynamic data is integrated with - the 
computer model. A feedback control scheme for the multiple- 
input, multiple-output system is determined using optimal 
control techniques. Time response results are obtained and 
analyzed. 

As a separate issue, the body of the device is modelled 
for electromagnetic analysis and a basic antenna design is 


determined for UHF transmission. 
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I. INTRODUCTION 


A. WHAT IS AROD? 

The U. S. Marine Corps has established the Ground/Air 
Tele Robotic Systems (GATERS) program to develop, test_and 
field robotic devices for direct use by U. S. Marine 
infantrymen. These robotic systems will be operated under 
ehe concept of tele robotics which implies that a man, 
stationed remotely from the vehicle, will control the device 
as if he is actually on board. 

The Airborne Remotely Operated Device (AROD) is the 
vehicle being developed through the GATERS program to meet 
the Marine Corps' needs for an airborne, short-range, direct 
support remotely piloted vehicle (RPV). The GATERS program 
office [1] states that, from a tactical perspective: 

The AROD is ... designed to allow the front line commander 
to see "over the next hill", out to a distance of two 
kilometers and "around the next corner" in an urban 
environment... 
This tactical requirement has led to a more technical 
statement of what AROD will be. The GATERS program has 
stated [2]: 
The AROD will be a ducted-fan, hovering device capable of 
a forward translation speed of 30 miles per hour. The 
vehicle will carry a fibre optic data link and on board 
cameras to support assigned missions out to a distance of 


5 kilometers. Normal operating radius is 2 kilometers. 


More detailed requirements are listed in Appendix A. 


do 


B. AROD AND OPTIMAL CONTROL 

The AROD is to be operated by a Marine, untrained and 
unskilled in aerodynamic flight. Therefore, AROD must be 
inherently stable and responsive to simple controls such as 
"Go up, dew; "rotate right/left", "O forward", "Cro 
sideways". These instructions will be commanded from a 
оу ск, through the control uplink iona ons a contio 
system. The control system will interpret the command and 
perform the necessary tasks to maneuver the AROD. The 


design of this controller is the major topic of this vemm 


Figure 1.1 is a schematic of the prototype AROD. 





Figure 1.1 Airborne Remotely Operated Device (AROD) 
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Chapter II discusses the AROD system and the equations 
which ~ its motion. The equations are complex and the 
system has multiple inputs and multiple outputs. One 
approach to controller design for such a — is to 
simplify the equations describing the system as: much as 


possible while retaining as much of the system characteris- 


“tics as possible. This is the method developed herein. 


The complex multiple-input, multiple-output AROD system 
is appropriately addressed using optimal control techniques 
with the aid of a digital computer [3]. Optimal control, 
therefore, is the approach taken and developed in this work 


to solve the problem of stability and control of the AROD. 


C. AROD AND RADIO FREQUENCY TRANSMISSIONS 

AROD, as previously described, utilizes a fiber optic 
tether to pass control information and a video signal from 
on board observation devices. Therefore, the need for radio 
frequency (RF) transmissions to or from the vehicle is not 
immediately obvious. In fact, under ideal tactical condi- 
tions,. the RF signature, susceptible to hostile electronic 
countermeasures, would be a most undesirable characteristic 
for a flying vehicle to have. Nevertheless, the fiber optic 
tether will be subject to damage and a back-up system must 
be considered for vehicle control, and for sending back the 
Video signal as well. Most importantly, training must be 


conducted prior to any conflict or situation in which AROD 
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is to be used. The fiber optic cable, not immediately 
reusable by the lower echelon ground forces, is expensive 
and, again, susceptible to damage. Therefore, it is 
desirable to provide a low cost method for training opera- 
tors and tactical commanders in the use of the AROD; hence 
the need for an alternative to a fiber optic tether. 

A back-up RF communications link to provide operator 
commands to the AROD has already been designed and imple- 
mented. However, this is an uplink only and provides no 
means for the on board cameras to return their wide 
bandwidth video signal to the ground station. Some 
consideration of this need has been made in the GATERS 
program and a transmission frequency in the 800 to 900 
megahertz range has been tentatively approved. [2] It has 
also been recognized that a video signal via RF might be 
adequate even though it was somewhat degraded relative to 
the high bandwidth video signals possible with fiber optics. 
With these — in mind, Chapter IV deals with the 
design of an on board antenna suitable for transmitting a 


video signal from AROD to the ground station. 


D. SCOPE 
In summary, this work will address two different 
problems associated with AROD: 
(1) Design of a control system for stable flight. 


(2) Design of a radio frequency antenna for the AROD. 
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1. Controller Design 
System design of the vehicle has been previously 
undertaken by AROD project engineers. The prototype AROD 
has been built and pertinent information is listed in. 
Appendix A. In Chapter II, the model describing AROD is 
developed. Design of the proposed control system is under- 


— 


taken in Chapter IIT. 
2. Antenna Design 
Chapter IV deals exclusively with design and 
analysis of an RF antenna proposed for AROD. The design 
approach is to analyze the requirements and determine a 
basic design which can be refined gradually until an 
acceptable solution is found. The computer program 
Numerical Electromagnetics Code (NEC), developed by the 
Naval besan Systems Center is used extensively to aid in the 
design. No attempt is made to recommend specific hardware 


components in order to implement the antennas. 
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II. AROD DYNAMIC SIMULATION 


The theory behind the equations governing aerodynamic 
flight is well developed. Procedures exist which make the 
modelling of complicated, multiple engine, highly maneuver- 
able aircraft almost a step by step task. Unfortunately, the 
dynamics which characterize AROD differ considerably from 
the standard flight model at a very basic level. Like all 
flying vehicles, AROD has six degrees of freedom of motion: 
three linear displacements and three angular rotations, but 
here the similarities end. This ducted fan device has no 
lift producing surfaces other than the propeller . which 
accounts for 100 percent of its lift and 100 percent of its 
complexity. This chapter traces the development of a 
suitable dynamic model of the AROD, concentrating first on 
the significant gyroscopic contribution of its propeller, 
and then on the complete vehicle. The goal is to providese 
computer simulation with which an experimental control 
algorithm may be integrated for the purpose of analysis and 


improvement. [4:pp.2-22] 


A. THE GYROSCOPIC AROD 
AROD is a gyroscope. The single propeller rotates about 


the longitudinal vehicle axis to produce a downwash or jet 
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of air through the duct which makes up the AROD body. This 
jet is directed by pitching the body forward, backward, or 
side to side so that the vehicle moves away from the jet in 
the direction desired. To appreciate the control effort 
required to pitch the vehicle in a specific direction 
without causing undesired effects (coupling) in the other 
five degrees of freedom, the equations governing a gyroscope 
need to be understood. 

Figure 2.1 pictures a gyroscope which rotates in the 
same plane as the AROD propeller relative to the coordinate 
system described. This axis orientation will also serve to 
define symbols for later use. For illustration, the body- 
fixed and earth-fixed axes do not coincide in Figure 2.1. 
X,Y,Z and P,Q,R are the body-fixed descriptive geometry for 
linear displacement and angular rotation rate. P,Q, and R 
will also be referred to as roll, pitch, and yaw. X',Y',2' 
and ao are the corresponding earth-fixed displacements 
and rotations. Х! сап ре thought of as altitude. 0,0, and 
y are the Euler angles. Velocities will be denoted by U,V, 
and W for rate of displacement in the X,Y, and Z directions, 
oem’ VW for X",Y',2”. Pictured also for use in later 
discussion are L,M,N (rolling, pitching, and yawing moments) 
СЕН си F (applied force components): 


Angular momentum, h, due to the propeller is defined as 


F e in, + jp + k. (2.1) 
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EARTH-FIXED 
AS 





BODY -FIXED hy, Fy, M 
AXIS 


Figure 2.1 Gyroscope and Axes Orientation 


The rotor lies in the YZ plane of the body coordinates. 
Since the coordinate axes have been chosen so that the rotor 
spins symmetrically about the X-axis, h, is directed only 


along X and h,,= h.,= ©. WEquallonw2 becomes 
В = Го =ih,, - . (2.2) 


At this point, a simple calculation will illustrate the 
significant influence the rotor has over the AROD. Refer- 
ence 4 states that in developing the equations of motion for 


aircraft with rotors which exert gyroscopic moments on the 


18 





body "more often than not, such gyroscopic moments turn out 
to be negligible." In the AROD case, however, 1, = .015 
ft*-lb(mass) and hover speed of the propeller is w = 7200 
rpm = 753.98 rad/sec. From Eq. 2.2, h,- 11.3 ft? -lb(m)/sec. 
Compared with a nominal total mass of 2.64 lb(m) (85 lb 
weight), it is clear that the angular momentum imparted by 
the rotor is significant and that gyroscopic effects will 
play a large part in modelling the re behavior of AROD. 

Complications arising from these gyro effects can also 
be illustrated. If the AROD body and rotor are considered a 
complete gyroscopic system, then the angular momentum of the 
body must be included in the formulation. Angular momentum 
Жоба косасіпа body is noted as hy — ihyy + Њу + ML 


With the orientation of Figure 2.1, Eq. 2.1 becomes 


h 


i (Myx, + hex) + jhi, + kh; (2.3) 


i (PI, + I,w,) % 101, % «RI, 


a mm 


where h is the combined angular momentum of the rotor 
amd body. 


Angular moments about the axes of a gyro are discussed in 


the next section and defined by Euler [5:p. 93] as. 


= h, = h,R E h, Q (2.4) 
M = h, - h,P + h,R 
N = h, - h,Q + h,P 


A simulation of this simple gyro model is given in 


Appendix B. Results are shown in Figure 2.2 for the case of 


19 


a disturbance about the Y axis (pitch). Note that the 


actual moments of inertia, listed in Appendix A, were used 


u 


RITON 25 
YAW -. 


ANGLE 
RATE 


(RAD/SEC) 





ГИНЕ 


Figure 2.2 Gyroscopic Coupling Effects 


in this first step computer ene The time response 
shown here is for the simple case with no rotor accelera- 
tion term (ш, = 0) and is the result of a one secon step 

for an angular moment, M, beginning at time = 1 second. The 
amplitude is unimportant and has no effect on the time 
response except to scale the result. The purpose of the 
input was to pitch over the coordinate system by applying T 


torque about the Y axis (Q). Ideally, a stable, well 


20 


behaved device would simply pitch over and maintain its new 
attitude. This is not the case for AROD and Figure 2.2 
illustrates the precession and nutation due to a severe 
coupling of the pitch(Q) and yaw(R) axes. These gyroscopic 
effects, then, are the reason a "simple" device such as AROD 
presents such a complex controls problem to the designer. 
B. EQUATIONS OF MOTION 

The basic equations of motion are developed for a 
typical aircraft in References 4 and 5. A combination of 
the two approaches is taken here to arrive at the set of 
equations describing AROD. 

The Euler equations of motion for a free flying rigid 
vehicle with six degrees of freedom are derived from 
Newton's second law. For angular moments applied to the 


coordinate axes, we have 


(2.5) 


e 

| 
a ja 
+Ë 


< 
|| 
a ja 
als 


N 


| 

= 
a ¡a 
+Ë 


where h= ih, + jh, + kh, and is the vehicle's total 
angular momentum. 


For linear or translational forces applied along the axes, 


_ ._dv _ dV Б ау 
F, = INO Fy = im +£ F, = kmt (2.6) 
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where v = iv + jv + kv and is the total linear transla- 
tional velocity of the center of mass(m) of AROD. 


The dynamic characteristics of AROD are represented by 
Equations 2.5 and 2.6. Other relationships are required to 
solve for each unknown and for each physical force and 
moment acting on the body. The development of a complete 
model is undertaken now for the specific case of AROD using 
actual measurements and — data from the proto- 
type. 

1. Angular Momentum 


Equation 2.4 is expanded to become 


h,= i(I,P+ Io.) + j (I, PR + I,w,R) (2%) 


k (I, QP + I,w,Q) 


-¡ I,RQ + jI,Q + kI,PQ 


D 
| 


h,= iI,QR - jI,PR + k1,R. 
where 
= _ on 
rei h =+ +06xh 
hz i (I.P + ТЕСИ ¡TO Е. 
w = iP + jQ + kR 


Combining Equations 2.4 and 2.7, the AROD angular moment 


equations are: 


Кој • 


L — I,P * (I, - I,)OR * 1,4, (2.8) 


O. 


M = I, 


I,R+ (I, - I,)PQ - QI w, 


+ (Ш - TI ET RO: 


N 


The rolling, pitching and yawing moments (L, M, and N) may 


be induced by the aerodynamics of the body and environment 
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еш ШШЕН и ры = through "the control vanes and are 
referred to as idile "applied moments". They are discussed in 
a later section. Equation 2.8 is next rearranged so that 
the unknown variables (dynamic states of EN rest) may be 


determined. Roll, pitch, and yaw rates are: 

BENE (хх E Y oR - ть (2.9 
STE [ ( y 220 г г ] | - 9) 
ee 

О бет ni o. waq 


R => І(І, - І,)РО + Ql,w, - N] 


Already obvious is the mathematical explanation for 
the coupling that was: observed in the last section. Every 
combination of pitch, roll, and yaw are present in Equation 
2.9. Pitch, roll, and yaw are each "coupled" to the others. 
Additionally, coupling of the rotor angular acceleration 
with roll will cause the throttle controlling the propeller 
speed to influence these states .and further complicate the 
system. 

2. Force Equations 
The expressions describing the translational forces 


are developed from the body accelerations which are 


v,= iU + jRU - «Qu (2.10) 
Vy= -iRV + jV + «PV 
V,- iQW - jPW + kW 


ere іу) + jv) + ky, and is the total body velocity 
weetor. 
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Equations 2.6 and 2.10 are combined and rearranged to obtain 


the differential equations containing force terms, 


U 


(RV - QW) + = | (2.11) 


У = (PW - RU) + ЗУ 


T 
| 


(QU. — By) + 


B 
N 


where m is the mass of the vehicle. 


Equations 2.11 contain nonlinear terms as did 
Equations sos These nonlinearities are in the form of 
products of states. The force components F,,F,,F, consist 
of lift from: the propeller (F,), other aerodynamic influ- 
ences (F,), and gravity TNT F, and F, depend largely on 
experimental data in this development and are taken up 
later. 

__ Gravity is a constant acceleration causing а 
constant force in the -X' direction (earth-fixed). Note 
that small changes in weight due to fuel loss are neglected. 
As the body rotates, the body-fixed axes change relative to 
the earth-fixed axes. Therefore, a coordinate translation 
is needed to relate gravity to each of the X,Y,Z amS 


(ЯАх,4у,42). The gravity components are given by 
i'g= ig, + ig, + kg, (2.12) 


Any translation from earth-fixed to body-fixed axes can be 
described by a rotation first about the Euler angle p, then 
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e, and lastly у. Figure 2.3 illustrates “these three 
rotations and the resulting relationships between coordinate 
systems. From Figure 2.3 amd Eq. 2.12, the gravity com- 


ponents are 


g, = g(cosy cosd) 2213) 
gy — -g(siny cose) 
: g, = g(sine) | 


The first-order differential equation for translational 


velocities can now be written as 


Ù = (RV - QW) +g, + E (2.14) 
: E 
V = (PW - RU) + g, + TY 
- Е 
W = (QU — PV) + g, ee 


We now have equations describing the six degrees of 
freedom of AROD and the foundation for the model has been 
vand. 

3. Kinematic Equations of Motion 

The coordinate translation through the Euler angles 
р,Ө,% resulted in an expression for relating gravity to the 
rotating body. However, the price paid for this necessary 
convenience is the addition of three more variables (p,0,+) 
to the equations. Therefore, three more equations are 
needed and, if possible, the three variables o,0,) should be 
expressed in terms of the existing force or moment vari- 
ables. Fortunately a relationship can be found in the 
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Earth to Body-F 


Figure 2.3 


equation describing the body angular rotations. Angular 
movement of AROD has been expressed in terms of the body- 
fixed angular rates P,Q,R. It can also be described in 


terms of the Euler rates o, Ө, апа y as 


о = ЛР + Ј)0О + кв = :'0) * j'O 9 k'y (2.15) 


“A coordinate transformation similar to’ the one 
illustrated in Figure 2.3 is possible to relate pitch, roll, 
and yaw with the Euler angle rates. However, some simpli- 
fications can be made with each rotation since o,0, are, in 
fact, the variables of interest. The resulting transforma- 


tion is called the kinematic equations of motion: 


Ф 1 sino cosptano P (2.16) 
Ө = О cosy -sinp Q 

а sino cose | 

у ` sino coso à 


— — 


4. Summary of Equations of Motion 

Nine equations and nine unknowns have been developed 
that completely describe the dynamic behavior of AROD. 
P,Q,R, and U,V,W represent the six degrees of freedom and 
9,0,», relate the constant force exerted by gravity to the 
body-fixed coordinates. Through the Euler angles, the 
behavior of AROD can also be observed relative to the Earth. 
Specifically, the translational velocities for the earth- 
fixed system U',V',W' can be determined from body-fixed 
velocities U,V,W and o,0,) using the coordinate transforma- 
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tion of Equatřiðn 2.13: The applied forces and moments are 
made up of external forces and commanded inputs, and are 
represented by sr, , wu and L, M, N. AROD can be 
pictured as the block diagram in Figure 2.4. The external 


influences are taken up next. 


C. APPLIED FORCES AND MOMENTS 

The applied forces and moments are a result of the 
commanded inputs and the movement of AROD through the air. 
The commanded inputs can change the rotor speed causing a 
change in the lift generated and a torque on the body. 
Commands can also come in the form of control surface 
displacements which steer the AROD and cause the body to 
rotate about its axes. The movement of AROD through the air 
will create drag. Drag forces can also apply both transla- 
tional forces and angular moments on the body and are a 


result of AROD's aerodynamic characteristics. 
1. Aerodynamic Characteristics 

Two methods exist for determining the aerodynamic 
forces and moments on a vehicle [4:p.111]: (1) by experiment 
Or empirical analysis and (2) by theoretical development. 
The great advantage of (1) 1S accuracy and the advantage of 
(2) is low cost. The fact that AROD is a ducted fan causes 
a complex interaction of the rotor with the body, and the 
downwash with the control surfaces. Hence, the ducted fan 


compounds the accuracy problem and empirical means for 
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determining the forces and moments is the most suitable 
approach. Sandia National Laboratories has built a proto- 
type AROD and has collected data based on wind tunnel tests 
of the AROD [6]. This data consists of tabular results 
describing the aerodynamic drag coefficients and downwash 
Swirl effect; stability derivatives; and physical measure- 
ments of constants such as weight, moments of inertia, and 
servo gains. This data is listed in Appendix A. The 
tabular data forms a basis from which applied forces and 
moments may be determined once the following flight informa- 
tion is known: 

(1) Weight ratio of actual AROD to 85 lb prototype. The 
data was originally obtained with an 85 lb prototype 
and must be scaled for other weights such as the 
model discussed here which is 76.5 lb. 

(2) Propeller angular speed. Propeller speed is 
determined by setting a hover speed (this model 7200 


rpm) and calculating changes from the throttle 
setting of the engine. 

(3) Propeller angular acceleration. Acceleration is 
based on the -difference between speed calculated 
from the throttle and speed at hover. This is an 
experimental approximation arrived at by  AROD 
project engineers. 


(4) Translational speed U,V,W. Aerodynamic forces on 
the body are based on the speed and direction of 
AROD (angle of attack). 


The forces and moments which result are of two 


types: aerodynamic and thrust. The forces and moments due 
to the aerodynamic data are referred to as Fi,, FF; VENE 
arc mi и ПЕ Those due to thrust and propeller con- 


siderations are Fy, Piy: Erz and L, MMe 
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No attempt is made in this work to verify the data. 
The relationships leading to the forces and moments dis- 
E in this section are taken verbatim from the informa- 

x _tion supplied by AROD project engineers. | 

| a. Angle of Attack 

x The tabular data (aerodynamic drag coefficients 
x and downwash swirl) result in forces and moments relative to 
x the total velocity vector, У,,,, of the body. A relation- 
ship is needed to transform these forces and moments to 
body-fixed axes forces and moments. Angles of attack (a, ) 
are the standard relationship for this purpose. Angle of 
attack (a) 1s normally defined as the angle in the XZ plane 
between V,.: and the X axis. Angle of attack (5) is the 


angle in the XY plane between V and the X axis (sometimes 


tot 
referred to as the side slip angle).[4] The angle of attack 
used in this model is defined and illustrated in Figure 2.5. 


The resulting relationships are given by 








sino = E (2597) 
тої 
š V 
Sing = 
Vet 
where Уон м ОМ + и 


b. Aerodynamic Forces 
The forces associated with the tabulated data 
ame lift (Fi), drag (F4), and side force (Е.).. They are 


also depicted in Figure 2.5. The transformation from lift, 
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Figure 2.5 Angle of Attack Illustration 


drag, and side force is accomplished in like manner to the 


transformation in Figure 2.3 and results in 


en 


„леве ! : 
Fs SlnoaCos8 -cosacosf sing Ei (2.18) 
г. = -sinasing cosasing cosa | [Е 
Fig -сова -sina Oo Ed 


C.  Aerodynamic Moments 
In similar fashion to the aerodynamic forces, 
the moments applied to the body axes as a result of AROD's 


movement through air can be derived. These moments are 
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referred to as the aerodynamic angular moments of roll, 


pitch, and yaw (L,, M,, N,) and are given as 


L, sinacosß -cosacosß sing y (2.19) 
M, = |-ѕіпоѕіпб cosasinß cosp В, 
| Fa -сова – 51па О Rp 

where 


Y, is the yaw moment relative to V.or, 
R. is the roll moment relative to Vot: 
P, is the pitch moment relative to V,,,. 
2. Forces and Moments for Control 

The forces and moments previously discussed are 
based on observation of an AROD prototype and are external 
to the system being designed. A second category is those 
forces and moments which are a direct result of the com- 
manded inputs. As mentioned uni itus. the inputs control (СІЗ 
the rotor speed and (2) the displacement of the control 
een. 

a. Moment Due to Ducted Fan Effects 

The AROD in hover (constant rotor speed) closely 

resembles a gyroscope. A gyroscope imparts no torque on its 
axis if it spins with a constant angular rotation. If the 
rotor accelerates (speeds up or slows down) a torque is 
applied to the axis. This torque is accounted for in 
Equation 2.8. However, AROD is also a ducted fan and the 


drag between the rotor tip and the inside body wall creates 


a moment about the X axis (P, roll). The project engineers 
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for AROD have determined an approximation for this moment 


based on experiment which is given as 
Lr. = -0729 (Aw.°) (ft-lbs) (2.20) 


where Aw, is the difference between the rotor speed and 
hover, Subscript Tr refers to a dependence on the 
throttle, and .0729 is a constant due to the duct and is 
referred to as K4,.. In Appendix A. 


— 


b. Moments Due to Control Surface Displacement 
Command inputs also include the ability to 


displace the four control vanes into the downwash from the 


duct. As these surfaces are displaced, they impart moments 
about the body axes. The vanes are arranged symmetrically 
as in Figure 2.6. The vanes are displaced by a servo 


mechanism connected directly to the top of each surface. 
Vanes (1) and (3) are operated together as "elevators" and 
impart a moment about the Y-axis Sim). Vanes (2) and (4) 
together are the "rudder" and contribute a moment about the 
Z-axis  (yaw). Vanes (1) and (3) displaced in opposite 
directions and (2) and (4) displaced oppositely work as 
"ailerons" to impart a moment about the X-axis (roll). 
Throughout this work, the term aileron, rudder, and elevator 
Will be used as defined here. The actual torque applied by 
each combination of vanes was determined experimentally and 
"Constants of effectiveness" were calculated by the AROD 


project engineers. 
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Figure 2.6 AROD Control Vanes 


The constants of effectiveness are referred to 
as dimensional stability derivatives when developed theoret- 
ically. They are given the symbols Le, Mee, Nre, for their 
contribution of moments about the roll, pitch, and yaw axes 
due to a displacement by the ailerons, elevators, and 
rudder. They can be thought of as acceleration quantities 
per radian of control surface displacement. The resulting 
relationships, concluding in moments about the three body 


axes, are 


35 


I, = йд (2.21) 


М, = М.б, 
N = Nreör 
where $ is the displacement of the aileron, elevator 


or rudder into the rotor downwash. Le, Mee, and N,e 
are given in Appendix A. 


Theoretically, each set of control surfaces will yield a 
constant of effectiveness about — axis (three constants 
for each aileron, elevator, and rudder). The symmetry of 
AROD causes negligible cross coupling of control surfaces, 
however, and it is ignored. Note also that the effective- 
ness of each control vane was determined for one rotor speed 
and that the change in downwash associated with the changes 
in rotor speeds is also ignored. 
3. Servo Equations for Control Surfaces and Throttle 

Commanded inputs sent to AROD include "go for- 
ward/backward", "go side to side", "turn around", and "go 
up/down". These commands can be satisfied by pitching the 
AROD forward or sideways (pitch and yaw) for translational 
flight, turning the AROD about its centerline (roll) for 
turning about, and changing the propeller speed for altitude 
control (throttle). Second order servo motors are used to 
displace each control vane as well as the throttle. 

a. Control Surface Servos 

The servos for each control surface are identi- 


cal. Three equations will describe the operation of the 
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elevators, rudders and ailerons. For each of these equa- 
tions, at least two servos are actually operating at the 
Same time on different control vanes. Each servo receives a 
command input (voltage) and results in an angular position 


for the servo. 


5, — -H46, - Hod, + Hou, (2522) 
Е $. — -HQ.6. - Hod. + HU, ü 
» = -H, 5. = Hs 6, zi H. u, 


where | 

H, and H, are servo gain constants, 

u, UE, and u, are the servo input. 

b. Throttle Servos 
The servo motor used to open and close the 

throttle is identical to the ones used for the control 
surfaces. An additional first order system describes the 
relationship between the throttle position and the power 
delivered to the rotor (thrust, T,). The complete third 
order system is 
бы = -H.5. - ње, + Hou, (2.28) 


T. = -T. + .5(5,) 


The power delivered to the rotor, T., directly relates to 
both the change in rotor speed, Aw,, and the force in the X 
direction. These relationships are 
ы. жет. е жғу, (2-24) 
АЕ. х= Т, 


where Pr,; 1S a constant relating the rotor speed to the 
power delivered by the engine. 
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D. ASSEMBLING THE MODEL: PROGRAMMING 

In this section, the equations of motion and the 
force/moment relationships are represented in a block 
diagram of the AROD model developed earlier, and the flow 


chart used to program the model is illustrated. 


1. The Complete Block Diagram 
" Figure 2.7 ils a block diagram of the complete model. 
For a more detailed explanation of the experimental aero- 
dynamic forces and moments, Appendices A and B should be 
consulted. 
2. Program Flow Chart 
The logic flow chart developed to program the 
na and the aerodynamic data is pictured in Figure 


Deer The IBM Dynamic Simulation Language (DSL/VS) source 


code used is given in Appendix B [7]. 


E. SIMULATION RESULTS 

Graphical illustration of time response is shown here. 
The four inputs u,,u,,U,, and u, are applied alternately as 
step inputs to the AROD model. A caution to the reader 
interpreting these results is appropriate here. The body- 
fixed angle rates P, Q, and R are the dynamic states of 
interest in the control problem. Emphasis should be placed 
on the response of these rates to control input and not on 
the earth-fixed Euler angles o, Ө, and y . The Euler angles 


are shown here only as a reference for the reader who wishes 
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Figure 2.8 AROD Simulation Logic Flow Chart 
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to picture what the vehicle's behavior would be relative to 


the Earth. Confusion can result if one tries to directly 
associate these angles with the roll, pitch, and yaw. By 
way of example (refer to Figure 2.1), the time derivatives 


Of о, O, and y are equal to roll(P), pitch(QO), and yaw(R) 
only when o - 0 and 9 and y are small. If the body has 
rotated 909 (о = 90), then Q actually approximates the 
derivative of y and R the СШ эуе of Ө. 
1. Zero Input Response 

Figure 2.9 is the resulting time response of the 
model when all control inputs are zero. For zero input, all 
control surfaces are parallel with the X-axis and not 
displaced into the downwash. The throttle is set for hover 
speed. Note that the pitch and yaw axes remain stable and 
zero as expected, but the roll axis spins opposite to the 
rotation of the propeller. | This precession is a result of 
downwash swirl which was measured and included in the 
tabular aerodynamic data. The propeller downwash is not a 
Miner flow of air but has a pattern of vortices which 
results in a force applied to the control surfaces leading 
to a moment about the roll axis. 

2. Aileron Input 

Figure 2.10(a) is the resulting time response when 
the aileron (6, ) is displaced and all other control inputs 
are zero as in Figure 2.9. Figure 2.10(a) indicates that 


the ailerons can be positioned to counteract the swirl 
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effect and actually roll the AROD oppositely. These results 
also affirm that the roll axis is uncoupled from the pitch 
and yaw axes as Equation 2.9 suggests. Figure 2.10(b) is an 
illustration intended to show that the roll — be  con- 
trolled and set to zero when no roll is desired. In this 
example, an initial displacement of 6,= -.0499 on the 
aileron control vanes is used and a bias equal to this 
initial displacement is added to =. 
3. Elevator and Rudder Input 

Figure 2.11 is the result when the rudder combina- 
tion of vanes (6.) is displaced into the downwash and the 
ailerons are set as in Figure 2.10(b) to minimize roll. 
Note the severe coupling between the pitch dti yaw axes, as 
expected. This coupling is a major complication to the 
control process and is dealt with in the next chapter. 
Figure 2.12 is the result when the elevator — И of 
vanes ($.) is displaced into the downwash, rudder is reset 
to zero, and the ailerons set as above. The results are 
similar to the rudder observations. In both cases, the roll 
has been affected by the pitch and yaw disturbances, a 
result of the tabular aerodynamic data which changes with o 
and f. The amplitude of the elevator, however, is greater 
than the illustration with the rudder and more clearly shows 


the serious problems that these aerodynamic moments present. 


44 





ROLL 
PITCH 
YAW 


EULER 
ANGLES 


4 








0 1 2 9 - 5 
ТІНЕ (ЕГ) 


u, *-.0499, u,=0, u,=0, ur=0 


Figure 2.10 (b) Time Response: Aileron Displaced 


to Eliminate Swirl Effect 


45 


ROLL 
PITCH 
YAW 


EULER 
ANGLES 


- 90 


u,=-.0499, u.=0, 


2 


3 


TIME (SEC) 


u, =.005, 





Figure 2.11 Time Response to Rudder Displacemnt 


46 


~ 


ROLL 


PITCH 
YAW 





ге нове 
© 
CLONE E. Щ 
------ У к 
“291 e. / 
2! ~ 
A 
C oog 
u === И 
Im „ЖЕ? 
а ~ à 
— PE. j CD 
| H 
= LN 
y “Р 
U) 
Qc uJ 
uJ — 
E. © 
D A 
ша 


(SEC) 
u, =0 


TIME 
u, *0,, 


-.0499, 1.105, 


a= 


u 


Time Response to Elevator Displacement 


Figure 2.12 
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4. Throttle Input 


Figure 2.13 is the result when the throttle (6,) is 
opened to increase the rotor speed and all other inputs are 
zero. The rotor applies a torque about the Fall axis and, 
indeed, this influence is seen in the results. The torque 
is directed oppositely to the spin direction (-P) and 
results in a less negative roll. However, the equations for 
pitch (Q) and yaw (R) include terms coupled to the propeller 
speed through the aerodynamic moments. The aerodynamic 
moments are determined by angle of attack which is a 
relationship between the body-fixed velocities. As U 
increases, the coupling with pitch and yaw becomes more 
pronounced, the angle of attack changes, and the moments 
determined from the tables change radically. The discon- 
 tinuity is best observed by the sharp change in roll, pite 


and yaw at .75 seconds in Figure 2.13. š 


5. Response Without Aerodynamic Moments 

It is insightful to observe the behavior of the 
model with the troublesome aerodynamic moments removed. The 
system behavior without the wind tunnel data supplying the 
moments is much more in keeping with expectations. Figure 
2.14 is the response of this modified system to a rudder 
displacement. Note the similarities of the nutation (Q and 
К) and precession (~) to the gyroscope illustration of 
Figure 2.2. Figure 2.15 is the response of the same system 


to an opening of engine throttle. This smooth system 
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response is clearly more in keeping with how we expect the 
actual system to behave. 

The validity of the data was not able to be con- 
firmed nor is it being challenged here. However, these 
results do suggest that the application of the data to the 
system model is a problem that is best addressed by the 


— 


engineers who derive it. 
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III. STABILITY AND CONTROL OF AROD 


The equations of motion describing AROD and the forces 
and moments which act on AROD were developed in detail in 
Gmepter II. A computer simulation of these ‘equations 
indicated solid stability in hover, but an inability to 
remain stable once disturbed from the hover state. There- 
fore, a method is needed to maintain stability and drive the 
system to a steady-state condition. Additionally, once a 
steady-state is achieved, it is necessary to steer the AROD 
up/down, forward/backward, side to side, and make it turn 
about its axis (roll). Hence, a method is also needed to 
control the AROD and cause it to change from a given initial 
steady-state to another desired steady-state. Stability and 
control are the topics of this chapter. Specifically, an 
algorithm based on optimal ir theory is developed in 
the form of the regulator problem to address the stability 
issue. The regulator is then adapted to the control issue 
formulated as. the tracking problem of control theory. 
First, however, a summary of the entire stability and 
control issue is offered. "em. a simplification of the 
complete AROD system is made so that standard, well founded 
methods may be used to address the stability and control 


dilemma. 
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A. SUMMARY OF THE CONTROL PROBLEM 

The task is to provide an algorithm on which to base a 
control system for AROD. The control system intended for 
AROD relies on inputs provided by an operator to command 
desired outputs. This section discusses the control problem 
based on the inputs, the outputs, and the equations which 


— 


relate them. 
1. Summary of the Complete AROD System 
The result of Chapter II was 18 first order dif- 
ferential equations which are listed in Table 3.1. Also 
listed are the equations for useful quantities which were 
derived in modelling the system. 
Figure 3.1 is a simple, open loop block diagram of 


the AROD system. The block labelled "AROD" represents all 





Figure 3.1  AROD Open Loop Model 


the equations listed in Table 3.1. The forces, including 
gravity, and moments acting on the block are external 
influences and are variables in the equations over which the 


designer has no control. The inputs of interest are u,, U, 
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TABLE 3.1 


AROD DIFFERENTIAL EQUATIONS 


| 
State 
È = Iy- OR- Io t+ h/t | 
Bee + O A N, + N, 1/1, < 
U = (RV - QW) + F,,/m+ F,,/m | 
V = (PW - RU) + Пау Шо Рука | 
W = (QU - PV) + F,,/m+ F,,/m 
ф = P + Osinp + Rcosptano | 
9 = Ocosp - Rsinp | 
) 7 Sino * Есозе 
| | 
Bla 
8 = -H,é, - Hoé, + Hou, 
$. = -Hi6. - Hod. + Hou, 
$ = -H,é, - Н6. + Ни. 
ôt = -Hy6, - Нб, + Но, 
u., апа ч, . These inputs are voltages applied to servos and 


determine the position of the control vanes in the rotor 
downwash and the throttle position. The outputs of interest 
are 9,9,%, и If these outputs can be manipulated at 
will, then the rest of the equations, which are derivatives 


of these, can be controlled as well. 
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2. Optimal Feedback Control 

Control systems may be divided into two basic 
categories: open loop and closed loop systems. Open loop 
control was attempted in Chapter II and is | in 
Figure 3.1. Step inputs for pitch(Q) and yaw(R) resulted in 
severe nutation and oscillations. Open loop control is not 
suited to the AROD application. Closed loop (feedback) 
control, then, is the direction taken in this work. 

A feedback control system is one which tends to 
Maintain a controlled variable at a fixed steady-state value 
[LONE Therefore, a control system must be able to sense 
variations in the System variables of interest in order to 
apply control to them. For AROD, the outputs m" the 
controlled system variables. A scheme to relate the output 
to the input is the feedback loop and is pictured in Figure 
5:02 The feedback matrix, K, must be determined so that a 
given control input always results in the same desired 
ОЦЕ Эа Note that the reference control Vector, =. 
not match the commanded input vector, Uu., except in steady- 
state. The time when r.# ц is called the transient 
response time and is also the time it takes the output to 
reach the desired state. It follows that transient response 
time should be as short as possible. Note also that the 
difference, r.-u., should be greatest when r, is first 
applied and should decay to zero as the outputs reach their 


desired state. This will result in a smooth transition of 
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Figure 3.2 Output Feedback Model 


the output states from beginning to end. In other words, 
there should be little or no overshoot of the desired state 
by the_transient states. A good feedback control will 
minimize both the transient response time and the overshoot 
of the desired state. 

A method of feedback control that might minimize 
both the transient response time and the overshoot is 
optimal control. Its name suggests that optimal control may 
provide the optimum solution to the problem of feedback 
control for AROD. However, it is more than semantics which 
drives the choice of controllers. 

First, optimal control lends itself nicely to a 


discrete time solution of the control problem. AROD will 


ov 


employ an on board, digital computer to perform inflight 
stability and control. While a continuous time controller 
can be easily discretized in many cases, design of a sampled 
data controller will simplify the procedure. | 

Second, optimal control is ideally suited to state 
variable techniques. A state-space approach is desired 
because of the enormity of the system and the complexity of 
equations. Matrix manipulation uiri ли и computers will 
prove far less tedious than traditional transfer function 
methods. 

Third, and by far the most important reason to 
select optimal control analysis over classical methods is 
the unique ability of optimal-control theory to accommodate 
multiple-input, multiple-output systems such as AROD [3]. 
Four inputs and four outputs have been identified. It will 
be shown in the next section that, although some separation 
of states is possible, severe coupling between other states 
precludes a reduction of the system to single-input, single- 
output systems. 

Fourth and last, optimal control is known to provide 
robust and insensitive solutions to the feedback control 
problem. Assuming that an appropriate performance measure 
is chosen to determine the optimal feedback gain matrix, K, 
the solution can be expected to have a fair degree of 
tolerance to plant model inaccuracies. This section 


develops an approximation of the full AROD model so that 
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troublesome nonlinearities can be dealt with. Therefore, 
some inaccuracies are bound to appear in the new, linear 
model compared with thes nonlinear system of Table 3.1. 
Clearly, robustness is not only a desired — of the 
controller, it is an absolute necessity if the controller is 


to be applicable to both the linear and nonlinear model of 


— 


ж. - 


AROD. 

With the advantages of optimal control comes a 
certain drawback: optimal control requires the availability 
ОЁ а11 system states [3:р. 22). Therefore, the feedback 
scheme pictured in Figure 3.2 must be modified to include 
all system states in the feedback loop. Figure 3.3 satis- 


fies this requirement. However, if the system states 


= = a y C) 





Figure 3.3 Full State Feedback Model 
are not available from sensor measurements or other means, 


computational observers must also be included to provide an 
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estimate of the hidden states. Observers add an extra 
computational task to the onboard computer, but otherwise 


present no problem to the control designer. 


В.  LINEARIZING THE AROD MODEL 

Throughout model development in Chapter II, nonlinear 
relationships were identified. These nonlinearities fall 
into two categories: nonlinear combination of states (e.g. 
Equation 2.8) and discontinuous functions (e.g. table look- 
up of aerodynamic force/moment coefficients). 

A linear approach to controller design is much preferred 
over a nonlinear analysis. Three major reasons drive this 
preference. First, the theoretical basis for controllers of 
linear systems is well developed and easier to implement. 
Second, in airplane design, the use of linear approximations 
for nonlinear systems has yielded good results for the case 
of steady-state flight and small disturbances therefrom [5]. 
Third, the optimal control solution, K, can often be reduced 
to a single, constant gain schedule for linear, time- 
invariant systems. A ‘constant gain matrix, K, will mean 
that a minimum amount of memory storage will be requis of 
the on board computer for control implementation [9]. 

In this section, the equations of Table 3.1 will be 
replaced with linear approximations based on the small- 
disturbance theory, steady-state assumptions, and selected 


physical approximations where necessary. 
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l. The Small-Disturbance Theory 


A straight forward discussion of the small distur- 
bance theory is contained in [5:p. 106]. Essentially, lf 
AROD is assumed to be in a certain steady-state reference 
condition, then motion of the vehicle consists of small 
deviations from the reference. If this steady-state is 
hover (all translational and angular movement is Zero) then 
each of the states in Table 3.1 will be zero or nearly zero. 
This is not.a bad assumption in the case of AROD. The 
flight condition in which AROD will spend most of its time 
is the hover. Only small deviations from the hover will be 
necessary to effect translational or angular movement. For 
example, a translational velocity forward (along the 2' 
axis) of four miles per hour requires only a 15 degree pitch 
of the body. The small-disturbance theory and the steady- 


state hover assumption result in the following simplifica- 


e- me 


tions: 
(1) Products of states are zero since the product of two 
small values is an extremely small value. 
(2) Derivatives of states are zero since the time 


derivative of a small valued time function is 
extremely small. 


(9) The sine of a state is equal to the state and the 


cosine of a state is equal to one. This is the 
small angle approximation for angles less than 15 
degrees. 


Table 3.2 is the set of equations which result when the 
above assumptions are applied to Table 3.1. Note that much 


of the coupling between states and all of the nonlinear 
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TABLE 3.2 


SIMPLIFIED EQUATIONS (SMALL-DISTURBANCE THEORY) 


State 
р =Р 
ONME 
- жан » 

ен ИЕ 
Q = -(I,/I,)u,R + (M + M)/I, 
R = (I,/I,)o,Q + (N, + N,)/I, 
о = с + ви 
V - F,,/m * F,,/m 
W = F,,/m + F¿¿/m 
$, = -H,é, - Hod, + Hou, 
6е = —Hi8, = Н + Ночь 
5. = -H,é. - Hj6, + H;u, 

Bii ó,  — -H46, - H6, + Huy, 


products of states have been eliminated. Note that the 
derivative term of rotor velocity, w,, remains in the 
equation for roll (P). The term is kept because ó, is not a 
State and because it is easily approximated. We wish to 
keep as much similarity between the linear and nonlinear 
models as possible, so physical approximations of 


troublesome quantities are important. 
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2. Physical Approximations of Forces and Moments 


‘The force and moment terms in Table 3.2 are made up, 
in part, of functions of aerodynamic coefficients measured 
at discrete AROD flight attitudes in a wind tunnel. The 
data, listed in Appendix A, follow no linear pattern and, 
instead, are discontinuous functions of angle of attack. 
The relationships used to include these data in the non- 
linear model were discussed in — II. If a steady- 
state condition is assumed, then some simplification of 
these forces and moments is possible. 

It is not desirable to exclude the aerodynamic data 
completely from the simplified model. An analysis of the 
magnitude of the forces and moments which result reveals 
that they are significant in all flight attitudes. However, 
if a constant hover is assumed then only one set of coeffi- 
cients need be considered. Hover condition is defined by 
AROD project engineers as о-909 апа 2-09 [6]. In hover 
condition, the forces and moments become those listed in 
Tapie 33. Note that, for the purpose of Table 3.3, the 
total body-fixed velocity is reduced to U. V and W are 
assumed zero. | 

Table 3.3 indicates that each force and moment in 
hover flight can be reduced to a linear function of a state, 
with one exception. Recall from Chapter II, we concluded 
that the rotor downwash created a swirl resulting ina ^ die 


on the control vanes and a moment about the roll 
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TABLE 3.3 


HOVER STATE FORCES AND MOMENTS 


Forces Moments 
Aerodynamic Aerodynamic 
= Т НО ТООС, 
- T.(.045) - 1.62 
Fagen U (бәс M. = Ue Sl) Q 
E. = =U (08) N,= -U(.18) 
Thrust Thrust 
F, = (76.5) + T.(2.125) L.- TP @5) — 6. (aes 
к,,= O М, = -5, (14.51) 
LE NO №. = -5, (16.68) 


Useful Approximations 
we> 753.98 rad/sec 


A (а алар 


axis. This moment was illustrated in Figure 2.10(a). This 
same constant moment appears in the equation for roll 
moment, L,, in Table 3.3. The nei factor of .0499, 
determined experimentally in the last chapter, can be 
derived by dividing the constant moment, --1.06, by the 
aileron control vane effectiveness, I4, (-21.29). This 
constant moment will be eliminated in the linear model, but 
Will be included in the resulting controller as a constant 
addition to the ailerem comtrol, wu The aerodynamic 


moments about the pitch and yaw axes are zero in hover 


flight. While this poses no problem to the linearizing 
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goal, it means that moments in other flight attitudes have 
not been approximated. Thus the controller resulting from 
the linear model will not consider these external moments. 
3. The Linear Model Assembled 
The complete, linearized equations are listed in 
Table 3.4. Numerical values for the constant quantities 


= 


TABLE 3.4 


LINEARIZED HOVER STATE EQUATIONS 


State 
| 
p =P 
2 =Q 
у =В 
P ="9(0296) -Т (3.738) - 1 (21.29) 
О = U(.4831) - R(6.75) - (14.51) 
R = U(.108) + Q(6.78) - $, (16.68) 
| Ü - T, (.894) 
ба  — -H.6, - H26, + Ни, 
W - “Hy 6, - Hod, + Hou, 


65 = -H,é, - Hod, + Hou 


EE 


ALL Q0 21.1. H. u, 


were taken from Appendix A. The next section will assemble 


the equations in a form suitable to determine the optimal 


gain matrix, K. 
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C. FORMULATING THE CONTROL PROBLEM 
1. State Space Representation 
The control system of Figure 3.3 is called a 
regulator. As discussed earlier, the goal of T€ regulator 
system is to obtain a stable closed loop system which will 
force any arbitrary initial states to zero. The first task, 


then, is to assemble the state equations in a form suitable 


for control analysis. The form desired is given by 
x(t)= Ax(t) + Bu, (t) (3 
where 
x(t) is the state vector, 
A is the plant matrix, 
B is the control matrix, 


u. (tc) is the control vector. 


Note that the A and B matrices do not vary with time. A 
time-invariant system greatly simplifies controller design. 

The 18 states listed in Table 3.4 are the complete 
states of the equations of motion for AROD. However, some 
of the states are unnecessary to the control issue and some 
additional assumptions will aid in determining the best 
control scheme. 

The control problem can be simplified if it is 
recognized Man translational velocities over the ground (V' 
and W'. in the earth-fixed horizontal plane) are accomplished 
by pitching the AROD forward/backward or side to и 
Therefore, the states important to this task are roll(P), 


pitch(Q), and yaw(R) along with their associated angles. 
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mer only translational velocity of interest is the earth- 
To vertical velocity or altitude rate. For small pitch 
and yaw angles, then, body-fixed V and W are equivalent to 
earth-fixed V' and W' and will be eliminated from the 
equations. Also, the earth-fixed vertical velocity (al- 
titude rate) will be approximated by the body-fixed veloci- 
22, ^- i = 
A caution about the — о,Ө,), is appropriate 
here. The linearized equations conveniently eliminated any 
dependence on the earth-fixed Euler angles. While this may 
be a valid assumption for determining the control solution, 
it changes the meaning of the angle states © and y. ф 
ылла» only on the angles © and ф which will be kept small 
in the hover condition, but © and y both depend on the 
Жай (һегеіге) angle, q, which can range from O to 3609. 
Therefore, the linearized angles ©, and yQ will be used T€ 
here on to remind us that these are body-fixed pitch and yaw 
angles and no longer relate the vehicle to the earth-fixed 
coordinate system. 
The complete, linearized system can now be put in 
the form of Eq. 3.1 and is listed in Table 3.5. Figure 3.4 
is the complete жоті flow graph for the linear system. 
2. Stability, Observability, and .Controllability 
a. Stability 
A gyroscope is a stable, open loop device. 


Small, bounded disturbances of the roll, pitch, or yaw axes 
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TABLE 3.5 


STATE-SPACE LINEARIZED SYSTEM 


x(t)= [0,0,;55,U;P) T, Riba rs. DO ПЕ EROS 


A(t) = 
ооо о 1 о о о O 0 0 0 0 0 0 0 
ооо о 0 9 Т 202778 0 0 0 0 0 0 0 
0 0: Өө Of OF (нс 0 0 0 0 0 0 0 0. 
ооо о о 28940 о O О 0 0 0 0 0 0 
O O O .024 O -3.74 O 0O -21.29 O 0 0 0 0 0 0 
000 07-0 -r 0 250 0 O 55 0 0 0 0 
000.4830 0 0 -6.75 0 -14.51 O O О 0 0 0 
000.1080 0 6.780 O O -16.68 0 0 0 0 0 
0000 о о ос 0 O О О О 1 0 0 0 
ооо о о о 0 9 O 0 0 0 0 j 0 0 
000 0 0 O 0 O 8 0 0 0 0 0 1 0 
ооо ооо о о OR @ 15 0 0 0 0 0 0 1 
ооо о о о о о -157.9 0 О O -17.77 0 0 0 
000 0 0 0 00 0O -157.9 0 0 0 -17.77 0 0 
850 60 0 о 0O QO 0 -157.9 0 0 0-17.77 Өң 
Роо0 o o o o 0 0 0 -157.9 0 O O -17.77| 
O 0-000- -0-0700 O 001157.9 0 0 0 је 
Bit) = 000000000000 0 157.9 O 0 
000000000000 O0 о 15.9 0 
оооооооооооо о O О 157.9 


Note: t is the transpose operator. 


produce bounded outputs. The outputs may be " the form of 
nutation ов oscillations ous pus will not grow with time . 
once the disturbance is removed. This stability was illus- 
trated for AROD in Figure 2.14 where the yaw axis was 
disturbed. The results compared favorably to the response 
of a gyroscope. 

The open loop linearized system may be tested 


for stability also. Since it is linear time-invariant, the 


68 


с нани NM MR oe TENER 





-0 
SIS 
N һо 
ro N 
Y. 
ds 
TN o 
“> 
Ф 
= 
Ф 
< 
d CH . a 
E c po 
а Y ” 
: ok e “O 
ЗЕ zu Tm 


6.78 


NS 


Figure 3.4 Signal Flow Graph of the Linear System 
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simplified system of Table 3.5 is stable if and only if the 
eigenvalues of the continuous time plant matrix, A, have 
real parts all less than zero [10]. Our system fails this 
test for stability since it results in У out of 16 
eigenvalues equal to zero. Therefore, the linear model has 
lost some of the characteristics of the nonlinear model. 
However, the control algorithm that stabilizes the linear 
system and minimizes the time denn and overshoot may 
still be applicable to the AROD. 
b. Observability 

Another relationship to test the applicability 
of the linearized plant to AROD is observability. Àn 
observable system is one for which all states may be 
determined from observations of the output. This is an 
important concept since an optimal controller will Бе 
useless if any system states are not known. To determine 
the observability of our linearized plant, the outputs which 
are avallable from direct measurement must be known.  Refer- 
ence 2 indicates that three single axis rate gyros, one- 
vertical rate gyro, a magnetometer (magnetic heading) and a 
barometric altimeter will be available on AROD. Control 
vane position based on a history of the input is also 
assumed to be known. Table 3.6 is a listing of the states 
which will be measured and the device used for each. 

The output of the linear, time-invariant system 


is given by 
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y(t) = Hx(t) (3.2) 


where 
u^ is the output vector, 
is the measurement matrix 
dimensioned o x n, 
n is the number of states, 
O is the number of measurements. 
| TABLE 3.6 _ 
MEASURED STATES 
Measurement Device State Obtained | 
Magnetometer e (heading) | 
Rate Gyro USED OD x 
Control History(Servo) 6, Or. h | 


and from Table 3.6 


Hi 

| 
O O O @ O O O Ə F> 
оооооооос 
ооооосооос 
оооооосоно 
ооо 
ooooo000 0 
Ооо ом о ОО 
O OOP OP OPOOO 
OO0OPOOOOO 
OOPOOOOOO 
OPOOOOoOOoOo O 
нн эәооооосос«осос« 
oOoooo0 00 0 
O00000 QO QO Oo 
оосоооосо 
ooo 00000 O 


Observability of a system due to its measured 
states can be determined from the observability matrix, W,, 
given by [11] 
Шь- EMEN , HAE . а... НА | (3.3) 
where 
ы ІШ айпшх mp matrix, 
m is the number of states measured. 


The rank of W, (the number of independent columns) must 


equal n if all of the states are observable from output 
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measurements. For the 16 state, linearized AROD system of 
Table 3.5, W, has rank 10. Although a control solution 
might still be found assuming all states were available, it 
would not be a realizable controller. Therefore, a dif- 
ferent approach must be taken. 

Figure 3.4 reveals that the roll(P) and the 
throttle (T, ) coupled together and the pitch(Q) and yaw (R) 
states are likewise coupled as expected from the gyroscope 
analogy. However, these two subsystems are coupled only 
through a force/moment relationship dependent on the total 
velocity of the body, M In the last section, V.,,. was 
approximated by U. A further simplification assumes that 
the moment on the pitch and yaw axis due to vertical speed, 
U, is negligible. The complete 16 state linear system can 
now be considered two independent subsystems of eight states 
each as pictured in Figure 3.5. State space representation 
Of the two subsystems is listed in Table 3.7. 

Each subsystem has an observability matrix of 
rank eight. Therefore, the states required in the subsystem 
configuration can be obtained from the measurements. 

c. Controllability 

A final test of the linearized model, before any 
effort is expended in actual controller design, is to deter- 
mine whether it is even possible to find a control sequence 
such that a desired final state may be reached. This con- 


cept is known as controllability [10:p. 205]. The control- 
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Figure 3.5 Linear Subsystem Signal Flow Graphs 
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TABLE 3.7 


STATE-SPACE LINEAR SUBSYSTEMS 


Roll/Throttle 


x (t)= [p Ar = [0 0.1 0 о 0 ТС 
U оо 0.894 0 о о о 
Р 0.020 -3.7 -21.290 о O 
| T. обоя о м оо 3 
| 6, 008p o 0 СР ШШШ: 
5, 0 070 o 0 Mo 91 
š, 000 0-157.9 0 -17.77 0 
5, ооо о 0-157.9 0-17.77 
Br =[o о о о о о 157.910 |5 
000 0 о о 0157.91 
Pitch/Yaw 
Xor (t)= [Oy Año = [oz o 0 O 
De оо о 1 O 0 0 0 
Q 00 0 -6.75 -14.510 0 0 
R 006.780 о -16.68 0 0 
5. оо о 0 O op 1 ип 
ó, ооо о 0 0 0 1 
бе 00 O 0-157.9 0-17.77 0 
š. оо о 0 0 -157.9 0-17.77 
Е Be = oto 0 5 O 157.9 O Tt 
00 о 0 O 0 157.9 
lability matrix is given by 
<= [В, АВ, А“В, ... ‚ А" В] (3.4) 


where 
W. is a nxin matrix, 
1 is the number of control inputs. 
The measure of controllability is the rank of W.. For each 


subsystem, the rank is eight. Hence, a control sequence can 


be found for the linearized system. 
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3. Discretizing the System 

a. Sampling Frequency 
We wish to obtain a discrete time, digital 
control algorithm. Therefore, the system of Table 3.7 must 
be converted to a sampled system. A basis for determining 
the sampling interval, AT, is needed. The sampling theoren, 


— 


based on the Nyquist criterion requires that [11]: 


1 
Ё. = — = 21, (3-5) 
where 

f, is the sampling frequency, 

AT is the sampling period, 

f, is the natural frequency of the systen. 

Multiple-input, multiple-output systems do not 

have a single natural frequency. However, if the system is 


considered a set of several single-input, single-output 
systems then a number of different bandwidths are observed. 
The highest bandwidth can be used as the natural frequency 
of the entire system for sampling purposes. In this way, 
four independent systems, with coupled terms set to zero, 
are extracted. The frequency domain transfer functions are 
listed in Table 3.8 along with their bandwidths. A Bode 
analysis yields the frequencies listed. Each of: the 
transfer functions also has the describing equation of the 
second order servo in the denominator. Second order 
analysis yields a servo natural frequency of 2 Hz. Since 


the servo frequency is the highest, it will be considered 
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TABLE 3.8 


SIMPLIFIED TRANSFER FUNCTIONS AND BANDWIDTHS 
Natural 
Frequency 


La e H2 
вё (S° +H, S+H; ) 


· 74 Ме 


Ме е Нг 


Б 5° (32 +H, s+H) ) -61 hz - 


Nr e H2 


= 5? (s?+H, s+H,) er 


2 299m 


= 3 (S* +H, StH, ) (St1) EL 





the natural frequency, f,, of the multiple system. The 
sampling frequency is now given as 


f, » 2f, - 4 HZ (3.6) 


Equation 3.6 is a lower bound on the sampling 
frequency that will allow the continuous time system to be 
reconstructed once sampled. Selection of the actual 
sampling rate is based on engineering judgment which calls 
for a factor of 10 between the highest frequency component 
and the sampling frequency [12]. A rate of 25 Hz meets this 
criteria and yields a sampling period of AT=.04 seconds. 

b. Analog to Digital Conversion | 

The process of conversion E illustrated in 
Figure 3.6. The block labelled D-A is a digital to analog 
signal converter. The block labelled A-D is an analog to 


digital signal converter. The AROD block is the nonlinear 
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Figure 3.6 Digital System Feedback Scheme 


model which has been linearized here. The next step is to 
combine the D-A, linear AROD model, and the A-D to obtain a 
discrete time representation of the linear, time-invariant 
system. Since the present linear model consists of two 
subsystems, the task is to convert the two systems and 
determine two control systems separately. The discrete 


version of the closed loop system is pictured in Figure 3.7. 


The discrete form of the system Equation 3.1 is given by 


x(k+1)= x(k) + Tu. (k) (3.7) 


where 
| k = KAT, 
k+1= КАТ+АТ, 
AT is the sampling period, 
15 discretized A plant matrix, 
P is discretized B control distribution. 
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x (K+) =ФхюЮ+ Го ® 





Figure 3.7 Discretized AROD Feedback 


The discrete plant and control distribution matrices are 


given by [11] 


Q = е^^Т | (3.8) 
Г 


where 
e is the natural logarithm operator, 
5 15 the Laplace operator, 
ds is the derivative with respect to s. 
Equation 3.8 can be solved by using a series 
expansion of the matrix exponential and integrated using a 
numerical scheme on a digital computer. The program called 


ORCONV [13] is used to obtain the discrete time system 


matrices of Table 3.9. 


78 





TABLE 3.9 


DISCRETIZED SUBSYSTEMS 


фьт E Tp; = 
10  .0400 -.0030 -.0167 0 -.0002 O -.0003 0 
01 0 .0350 0 .0004 O 0 0 . 00001 
O .0009 1  -.1466 -.8217 -.0015 -.0134 0 -.0299 -.00003 
00 0 .9608 0 „0189 0 .0003 0 . 0007. 
BO 0 0 .9010 0 .0275 0 .0990 O 
00 0 0 0 .9010 0 .0275 0 .0990 
00 0 0 -4.334 0 .4132 0 4.334 0 
00 0 0 0 -4.334 O 4152 0 4.334 
Dor Е Це = 
1 0 .0395 -.0054 -.0113 .0012 -.0001 .00001 -.0002 .00001 
O 1 .0054 .0395 -.0010 -.0130 -.00001 -.0002 -.00001 -.0002 
0 0 .9636 -.2667 -.5530 .0879 -.0090 .0010 -.0202 .0016 
0 0 .2681 .9636 -.0768 -.6357 -.0009 -.0104 -.0014 -.0234 
оо 0 0 .9010 0 10275 0 .0990 0 
оо 0 0 0 .9010 0 .0275 0 . 0990 
00 0 0 -4.334 0 „4152 79 4.334 0 
O 0 0 -4.334 O . 4132 0 4.334 


D. THE OPTIMAL REGULATOR 
A review of the linear control problem is appropriate 
here. We now have two discrete-time subsystems, each. 


defined by Eq. 3.7. The roll/throttle subsystem is descri- 


bed by the plant matrix Dp 7 and” control distribution Гь.. 
The pitch/yaw sust is described by qu. апа Гор These 
matrices were listed in Table 3.9. We also know that all of 
the states in each subsystem are available either by direct 
measurement or computational observers. We wish to find the 


closed loop feedback matrix, K, in Figure 3.7 such that any 
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disturbance of any state from the reference hover condition 
(all states zero) will be damped out, returning the states 
to zero in minimum time with minimum overshoot. Since each 
subsystem is controllable, the control law that satisfies 


these requirements is given by [10:p. 771] 


u, "(k) = -Kx(k) EC 
where " denotes optimal in the minimum time, minimum 
overshoot sense. 

Note that K is a constant matrix. We state this goal to 
reduce to a minimum the amount of computer memory required 
to store the feedback gains and, in general, to simplify the 
control algorithm. The validity of this approach is born 
out by experience with linear time-invariant systems [3]. 
The applicability of the resulting controller to the 
nonlinear system, like the previous linearizing assumptions, 
Will be determined by the time response of the full AROD 
model once K is found. The following paragraphs will 
provide analytical justification for the steady-state 
approach. 
1. Optimal Control 

Optimal control theory is concerned with minimizing 
a us criterion chosen by the designer [3:p. 3]. 
аса the primary focus in developing the controller 
should be on a performance measure that quantifies the goals 
of the design effort. These goals are summarized here: 


(1) Minimize the transient response time 
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(2) Minimize the state overshoot 


(3) Determine a constant gain schedule, K. 
(4) Operate within the physical constraints of the 
system. . 


The performance measure typically used in formulating the 


discrete time optimal regulator control solution is 


cM 
И. 


= -X(N)' Cx (N) (3.10) 
N- 1 

а > Y [x(K)'O(k)x(k) -* u, (K)* R(k)u, (X) ] 
k=0 


J is the cost function, 

N is the final time step (NAT), 

k is the time step index, 

Q(k) is the state weighting matrix, 

Е(К) 15 the control cost weighting matrix, 
e is the final state weighting matrix, 
€ is the matrix transpose operator. 


The proper choice of C, Q(k), and R(k) will ensure that the 
= schedule, u,(k), that minimizes J will meet the 
füset р" of our design goals and minimize the control 
effort required. 
The solution to Eq. 3.10 for time-invariant systems 
is found from the recurrence relation 
K(N-k) = -[R(N-k) + TEP(k-ı)T] ! (3.11) 


х ГЕР(к-1)ф 
апа 


P(k)= M(N-k)*P(k-1)M(N-k) + K(N-k)*R(N-k)K(N-k) + Q(N-k) 


and 
м(ҝ) = ф - TK(N-k) 
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To meet goal (3), a constant K, it is necessary to 
know that under certain conditions the recurrence relation 
converges to a steady-state for large enough К. The 
conditions sufficient for a steady-state convergence are: 
constant, positive definite weighting matrices Q and R; 
final state weight C=0; and u. (k) unconstrained [10]. The 
final time loses its meaning for the steady-state case since 
we assume convergence prior to reaching time-step N. 
Equation 3.10 becomes 

Кс 
J= Y [x(K)'Ox(X) -* u, (K)* Ru, (X) ] (3.12) 
k-0 
where k, is the time step when J converges. 
The recurrence formula can now be written [13] 
к t -Tah 
K = [R + l' PT] 'T* PÓ (3.13) 
and 
Р = МРМ + КВК + НОН 
апа 
"a M = - TK 
where the measurement matrix, H, modifles Q so that only 
states of interest are considered. 

The revised performance measure appears to have no 
dependence on the final state. In fact, the solution, K, is 
sometimes referred to as a "suboptimal" gain schedule. 
However, experience shows that for time-invariant systems, 
the suboptimal gain schedule is often quite satisfactory. [9] 

2. Selecting the Performance Measure 


The weighting matrices Q and R must be selected to 


define the cost function used to determine the optimal 
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gains. For the linear AROD problem, Q and R are both 
diagonal matrices. Each diagonal element of Q corresponds 
to a weight placed on the deviation of a particular state 
from its desired value. Equation 3.13 introduced the idea 
of weighting only certain states utilizing the H matrix. 


The states we are interested in weighting are the angle 


= 


states, the altitude rate, and the control vane displace- 
ments and displacement rates. The angle states and altitude 
rates describe the steady-state of the system. The control 
vanes are constrained in their movement. Therefore, these 
are the states with which we are concerned with controlling 
through the gain schedule. A definition of the measurement 
matrix, H, strictly for the purpose of the cost function, is 


(3.14) 


ооооонм 
оооонро 
oooooo 
oooooo 
ооороо 
оонооо 
омоооо 
»"ооооо 


where H applies to both subsystems. The non-zero 


elements correspond to O > for-roll/ throttle 


ӚТЕ O, ¿Ue ro +82 ¿8% ‚s. for pitch/yaw. 


Now Q is a 6x6 matrix. If all the states of interest are 


weighted equally, Q is given by 


(3.15) 


Ooooonm 
оооосҥн+о 
ОоОооноо 
OQ OF OOO 
OPOoOoooO 
MOoOooo6o 
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This reduction in cost states will streamline the iterative 
design process of the gain matrix in the next section. 

Each subsystem has two control inputs: u, and u, 
for roll/throttle, and u., и, for pitch/yaw. A control 


weight matrix which treats each control cost equally is 


= |1 0 (3.16) 
М ы E 4 | 


The cost function resulting in the best gain matrix 
for our purposes may not be the one which weights each state 
and cost equally. Factors seh as the units a particular 
state is scaled to and system constraints on other states 
require scaling the diagonal values of the Q and R matrices 
accordingly. It may also be more important to drive a 
particular state more rigidly than others because of the 
dynamics of the system. In the case of AROD, it is most 
important to drive the angle states and the altitude rate to 
the desired state because the other states, all derivatives 
of these, follow. The time response of the angle states and 
altitude rate is the standard by which the best cost 
function is selected. Therefore, an iterative approach is 
taken, altering Q and R as necessary to meet the following 
specific objectives: settling time of two seconds to less 
than 10% overshoot, and maintaining constrained states 
within limitations. 

The settling time requirement is chosen for two 


reasons. First, the servo response time (second order time 
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constant) is .5 seconds which places a lower limit on the 
response. The equations are coupled which means action by 
one state will force a delayed, but proportional reaction by 
the others. Two seconds for any response on a slow speed 
hovering vehicle such as AROD is considered adequate. 
Second, if the erratic behavior observed in Chapter II is to 
be controlled, then a quick response to disturbances and 
commands is necessary. 

The control vanes are constrained in displacement to 
1.5236 radians (30°) and in velocity to +.8727 radians per 
second. The throttle is constrained to +100 radians and 
+100 radians per second. The iterative design process must 
account for these constraints by weighting the corresponding 
elements of the Q matrix. Appendix A lists constraints. 

Constraints on the servo states imply that the rate 
of change of other states may be restricted as well. A 
constraint on the output states means a limit must be 
applied to the disturbances from which we expect the system 
to recover in the given time. An angular displacement of .1 
radians (6°) is a IEEE ccs disturbance to the system. 
This does not mean that the system is restricted to motions 
within 69 of the hover condition, only that the response 
time for which we design applies up to these limits. 

4. Iterative Design Process 
We are now ready to form cost functions for each 


subsystem and determine suitable gain matrices. A numerical 
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algorithm for solving Equation 3.13 is discussed in [3:83]. 
A computer program for iterating the solution and converging 
to the steady state is described in [13] and is used to 
determine the gain matrices in this section. 

Once a gain matrix is determined, the time response 
of the linearized, closed loop system can be tested using 
the linear system computer simulation in Appendix C. The 
process repeats until an acceptable time response is 
obtained. 

а. Roll/Throttle SubSystem 

Table 3.10 summarizes the results with various 
cost functions. The values for Q are the diagonals of the Q 
matrix which are the weights on p, U, 6,, &, 5 қ апа бей 
respectively. Тһе values for R are the weights on u, and u, 
in that order. 

Iteration (1) is with a Q normalized to account 
for scaling differences between states. Results indicate 
that the normalized cost и is not adequate. A 
shorter response time is desired. 

Iterations (2), (3), (4), ¡and (5) chronicles 
trial and error design process where each succeeding cost 
function is chosen based on the response of the previous 
ones. A dominant consideration is the coupling between the 
throttle and roll states. Since is undesirable to have the 
throttle used to correct the roll angle, a heavy cost is put 


on maintaining the throttle at the desired state. 
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TABLE 3.10 


ROLL/THROTTLE SUBSYSTEM COST DESIGN ITERATIONS 
Settling Overshoot 


Cost Initial Time % (absolute) Deviation 


(1) =.1 ae (0) .02 0 .2 


(2) 14 1.5 (0) „03 -.0002 .35 


(3) š | (0) .03 -.004 .35 
(4) í (0) .07 -.0007 1.2 


(-.0001) .055 .0125 .87 
1. (.004) .045 -3. .275 


Cost Functions 


(3) 
R 


Q Q Q 
20 50 P 50 
60 1000 2000 11612000 


«OL ¿01 


3 2 
ker а 





|. The final iteration results in an acceptable 
time response and satisfies the state constraints. The 
optimal gain matrix, Kp; converged in 14 iterations at 
AT=.04 seconds рег iteration. Hence, the steady-state 
constant gain schedule is reached after .56 seconds. Since 
the settling time is 1 and 2 seconds, our assumption that a 
constant gain schedule would provide the optimal solution is 
valid. 
b. Pitch/Yaw SubSystem 


Table 3.11 lists results for the iterative 


design of the pitch/yaw subsystem gain matrix. The values 
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for Q are the diagonal elements corresponding to Өр, Yb; der 
ORG js апа 5, respectively. u, and u, are the control 


inputs in R. 


TABLE 3,11 


PITCH/YAW SUBSYSTEM COST DESIGN ITERATIONS 
| Settling Overshoot 
Cost Initial Time % (absolute) Deviation 
E Value Ө ф Oo  » “s OM: 


1) Ө 

(2) © 

(3) Ө =.1 .75 1.5 3% (.0125) .05 -.065 .7 -.7 
Y 


=.1 эһ 1.5 10% (.013) .035 -.055 .35 -.45 


1 абы. 1: 8% (.0125) .08 .045 .6 .5 


0 1. 1.5 (-.012) 2% .065 .06 .75 .6 


| Cost Functions 
(1) (2) 
Q R Q =R 
100 1 1 
100 1 aL 


„ Ой 





The same procedure as roll/throttle for iterat- 
ing through various possible cost functions is used for the 
pitch/yaw subsystem. | As in Table 3.10, only selected 
iterations are listed. Insight gained from the roll/ 
throttle iterations was valuable in reducing the steps 
needed for the pitch/yaw design. 

The final gain matrix, Kaos, converged in 13 
iterations, or .52 seconds which is less than the settling 
time of .75 seconds. Therefore, the steady-state gain 


matrix is an optimal solution for this cost function. 
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c. Steady-State Gains 
Table 3.12 lists the optimal steady-state gain 
matrices resulting from the iterative design. A test for 
stability places all eigenvalues of the айе time closed 
loop subsystems inside the unit circle. Therefore, the gain 
matrix, K, has succeeded in stabilizing the open loop plant. 
TABLE 3.12 


STEADY-STATE GAIN SCHEDULES 


Roll /Throttle 


-1.99 н .274 1.75 .00994 .140 .00016 


2385 š Done 25529 1.67 60962 .137 


.320 1.05 -.208 .127 -.00264 


-.266 ¿107 _T.E3 | ¿00096 .150 





5. Linear System Optimal Requlator Response 

The graphical results discussed here are for the 
optimal gains of Table 3.12 applied to the linear model with 
state constraints. The DSL/VS simulation is listed in 
Appendix C. 

Figures 3.8(a), (b), (c), and (d) are the regulator 
time responses for initial conditions of g(0)=.1, U(O)=.1, 
©, (0)=.1, and y4(0)-2.1, respectively. The design goals for 


overshoot and settling time are satisfied for the linear 
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Figure 3.8 (a) Linear System Regulator Time Response 
with Optimal Gains, p(0)=.1 


90 


PITCH 

YAW 

ANGLES > 

CIO? RAD) Tr оь 


ROLL 
ANGLE 
(16? ВАО) 


ALTITUDE 
RATE _ 
(5 XI FTAEC) 





2 3 
ТІЛЕ (SEC) 


Figure 3.8 (b) Linear System Regulator Time Response 
with Optimal Gains, U(0)=.1 


91 





PITCH 
YAW 
ANGLES 
(10° RAD) 


ROLL 
ANGLE 
CIO? RAD) 


ALTITUDE 
RATE. — 


(5 X10" FT/SEC) 


5( 


50 


Figure 3.8 (c) 


-д, 
NG 
NY 
sy 
1 2 3 5 
TIME (SEC) 
peu 
4 5 


2 3 
ТІНЕ (SEC) 


Linear System Regulator Time Response 
with Optimal Gains, Ө(0)-.1 


92 











122г5-------------------------------------------- 
ICH \ 
PITCH y 
YAW 4 \ b 
ANGLES | | 
(16? RAD) La 
HA O nn nn nn nn 
2 
O, 
“So 
157 Time (SEC) 
ROLL 
ANGLE jm O roces 
CIO" RAD) 
ALTITUDE 50 
RATE — 
(5XIO?FT/SEC) у | 
“Sa 
0 1 4 © 


e 3 
РМС (SEC) 


Figure 3.8 (d) Linear System Regulator Time Response 
with Optimal Gains, y(0)=.1 


95 


system. These results, for which the optimal gains were 
designed, are a reference to which results with the non- 


linear system can be compared. 


E. THE OPTIMAL TRACKER 

The optimal gains derived for the regulator can be 
extended to the tracking probleme тие tracking probie 
characterizes the command-response relationship desired to 
steer and maneuver AROD. In theory, several approaches to 
the tracking problem exist. The method used in this work 
utilizes the error state. Тһе control law which applies is 

и. (к) = -K(x(k) - r¿ (k)) (3.17) 


where г; (К) is the reference (commanded) state vector. 


This control law is pictured in Figure 3.9. The error state 





Figure 3.9 Tracker Control Scheme 
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as it drove the system states to Zero in the regulator 
problem. A full state reference vector is required which 
means a commanded value must be specified for every system 
state. The reference vector for the rollt subsystem 


becomes 


r - Гес,0с,Р,Т, sôa dt „ба „ „бе 1" (3.18) 


e. Срт с 


and for pitch/yaw 


r = [Oe Pe Qe Re re 6, „бе „ »бе „ Ј" (3.19) 


сор с 


where the values of the reference vectors may be time 
varying and are the desired values of the final state. 
Application of Eq. 3.17 to the linear model is described 
in Appendix C. Figure 3.10(a) is the response for desired 
states ф=Ө=.1. The response is within the desired limita- 
tions for settling time (2 seconds) and overshoot (10%). 
Coupling is quite pronounced between the pitch and yaw. 
Figure 3.10(b) is the response to U-.5 and coupling between 
the throttle and roll heading is seen. Finally, Figure 3.10 
(с) is the response to a step input on yaw (ф=.1). The 


results are within specification. 


F. RESULTS WITH THE NONLINEAR SYSTEM 

The optimal steady-state gains for the linearized AROD 
model are determined in the previous section. The goal in 
designing the control algorithm is to provide fast response 
to disturbances and command inputs, while maintaining 


minimum overshoot of the final state. The requirement for 
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quickness is driven by the erratic behavior of the full 
model in Chapter II when it is disturbed from the steady- 
state pitch or yaw positions. This section applies the 
linear system controller to the nonlinear —— Results 
are given first for the controller applied to the nonlinear 
model without the pitch and yaw aerodynamic moments and 
second to the full model with all moments and. forces. 
1. Nonlinear Model without ——— Moments 

Figure 3.11(a) is the nonlinear system step response 
for the heading angle, o. The settling time is 1.5 seconds 
which is the same as for the linear system. However, the 
overshoot is nearly 15% which is larger than the linear 
case. Figure 3.11(b) is the response to a desired altitude 
ree U! of .1. Note that for the nonlinear model, the 
altitude rate is derived from the Euler angles and it is no 
longer approximated by the body-fixed velocity, U. The 
altitude rate settles in less than 1.5 seconds, but the 
heading angle, disturbed through coupling, settles in 1.75 
seconds with .05 radians overshoot. Still, the system does 
settle quickly and with zero steady-state error. 

Figure 3.12 is the system résponse for a desired 
State step input on yaw, dy, of .1. Severe coupling to o 
and U' is observed, while pitch, 9,, behaves as expected. 
The reality of the linearizing assumptions becomes evident 


here. The results simply indicate that as the body is 


99 


PITCH 
YAW “° 30 
ANGLES —=— Oy Ye 
СЮ? RAD) Г 4 
-50 
0 1 2 3 4 5 
- TIME (SEC) 
ROLL m s. E 
ANGLE 
сю? RAD) -—® 
ec 
ALTITUDE | 
RATE | | U 


510? FEO e — — ————— 


2 3 
TIME (SEC) 


Initial States: 2=0, Џ'=0, 0,=0, ,=0 


Desired States: р=.1, U'=0, 0,=0, %,=0 


Figure 3.11 (a) Nonlinear System Tracker Response 
Heading Angle Command, o 


100 


PITCH 
YAW 
ANGLES 
CIC? RAD) 





e 3 
TIME (SEC) 


ROLL 
ANGLE / 
CIÓ? RAD) FU 


ALTITUDE - № m 
RATE , | 
(ХО РУЗЕС) _ || 


10 — — un 01 XL 


2 3 
TIME (SEC) 


Initial States: о-0, U'=0, &,=0, %p=0 


Desired States: р=о, U'-.1, 0,=0, ,=0 


Figure 3.11 (b) Nonlinear System Tracker Response 
Altitude Rate Command, U' 


101 








190 =a q A НИ 
PITCH / 
Í 
YAW 
- „| y 
ANGLES | b 
(ЮЗ ААО) 4 
ROLL еј И 
ANGLE — 
(ЮЭ ВАО) 
ALTITUDE 5 
RATE m 
БХК | 7 —. 2... 
Т 
CU 
50 \ 
0 1 4 5 


2 3 
TIME (SEC) 


Initial States: «(p=0, U'=0, 6,=0, ,=0 


Desired States: p=0, U'=0, Өр-0, фь=.1 


Figure 3.12 Nonlinear System Tracker Response 
Yaw Command 


102 


pitched (or yawed), some of the altitude holding thrust will 
be lost and redirected to a horizontal velocity. Figure 
3.13 shows that if the pitch or yaw is kept small enough 
(here ,=.03 radians), then the altitude rate can also be 
controlled. Figure 3.14 goes a step further and illustrates 
that if a bounded input (pitching the body for 1 second) is 
applied, then a speed across the ground can be realized 
while maintaining the desired heading and altitude rate. 
2. Nonlinear Model Complete 

Results for the controller applied to the full model 
including the troublesome pitch and yaw aerodynamic moments 
are shown here. Figure 3.15 is the system response to a 
step input on roll (p=.1) A steady-state error of -.04 is 
and is due to the downwash swirl moment previously dis- 
cussed. However, сина ЕО the other states causing the 
instability noted in Chapter II has been eliminated by the 
"uer. Figure 3.16 is the result of a step input for 
U'=.5. Coupling of the velocity to the pitch and yaw axes 
causes enough of a disturbance in angle of attack to upset 
the system beyond the constraints of the control vanes and 
the controller is unable to correct it. To illustrate the 
behavior of the data to pitching and yawing, Figure 3.17 


shows the resulting aerodynamic moments, L,, M,, N,, for 


step input for altitude rate. 
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3. Conclusion 

Clearly, the controller is not adequate to deal with 
such large discontinuities as the aerodynamic moments 
present. The control vanes are the Net ing factor. 
Therefore, the data or the way it is related to the states 

in the model is not correct. 
` However, we concluded in Chapter II that the 
behavior of AROD without these moments closely resembles the 
gyroscopic model that we expected. If the controller is 
judged based on its performance with the nonlinear model 
without the pitch and yaw moments, then the control design 
is valid. The fact that the controller maintains stability 
while directing commanded inputs to steer and maneuver the 
vehicle supports the linear model approach taken in this 


work. 
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IV. RADIO FREQUENCY ANTENNA DESIGN 


A. REQUIREMENTS FOR A RF ANTENNA 

Detailed technical requirements for the AROD video 
downlink have not been specified by the GATERS project 
office as of this writing. Neither have the tactical or 
operational requirements been addressed in detail. However, 
GATERS project quarterly review conferences have undertaken 
discussion of this aspect of AROD and some consensus has 
emerged [2]. It is from these discussions and from the 
stated operational mission of AROD (see Chapter I) that 
assumptions are made on which to base the antenna design. 

Several facts are known about AROD which pertain to RF 
wave propagation. 

(1) Cost and weight are key factors and a simple design 
utilizing inexpensive, readily available materials 
is most appropriate. 

(2) The AROD is intended to be man carried and easily 
readied for launch. Rugged, lightweight 
construction is required. 

(3) A frequency range in the UHF band of 800 to 900 
megahertz (MHz) has been tentatively identified for 
video signal transmission. 

(4) A range of at least five kilometers is required if 
AROD is to be operated at the maximum planned 
distance from the operator. 

(5) An omnidirectional pattern is desired for training 


purposes to ensure that the video signal can be 
received regardless of AROD's orientation. 
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(6) Angular speed of the propeller during hover is about 
7200 rpm (120 Hz) and a three blade, carbon fiber 
composite propeller is to be used. 


(7) The body is made of a carbon fiber composite and the 
frame is composed of aluminun. The forebody is 
fiberglass. 

(8) The AROD contains an on board computer, sensors, and 


a radio receiver for the command uplink. 
This list by no means includes all of the factors which 
contribute to the electromagnetic character of AROD nor does 
it address every possible requirement for a RF antenna on 
AROD. However, these requirements allow a design goal to be 
determined based on the following assumptions: 


(1) The 120 Hz disturbance from the rotating propeller 
and the command uplink RF signal will not interfere 


with the 800 to 900 MHz transmission signal. An 
order of magnitude between signals is a common rule 
for non-interfering signals. The propeller (7200 


rpm) and the command uplink (72 MHz) meet this 
constraint. 


(2) For the purpose of electromagnetic analysis, the 
AROD can be modelled as the cylinder pictured in 
-— Figure 4.1. The cylinder has the outside 
dimensions of the AROD main body surface shown in 
Figuuses 1«1. 


(3) The carbon fiber body is highly conductive. AROD 
project engineers have stated that the composite 
used has high conductivity characteristics. 


(4) The carbon fiber body will essentially shield out 
electromagnetic interference from the on board 
computer and sensors. 


(5) The antenna will conform to the body of AROD as much 
as possible to prevent breakage while being carried 
and snagging during operation in closed-in areas. 


(6) The antenna will be fixed to facilitate simple 
launch preparation. 
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Figure 4.1 AROD Model for Antenna Simulation 


B. DESIGN SELECTION 
A -togical point to start the choice of antenna is with 
devices already in use on similar vehicles. AROD is most 
Similar to a satellite in that it is free flying yet not 
fast почне and aerodynamic influences are negligible. One 
UHF satellite antenna in common use is the mesh array 
commonly referred to as the umbrella [I$ = However, the 
umbrella antenna is neither rugged nor omnidirectional and 
would be a poor choice for AROD. 
| Stripline antennas are arrays of slotted strips which 
are surface mounted to a vehicle. The slotted strips are 


cut for specific frequencies and are generally not tuneable. 
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These antennas are lightweight, able to be flush mounted and 
could be made rugged and durable for AROD. Research for 
this report has found only L-band and S-band variations of 
the stripline satellite antenna in use. However, the 
concept of a flush mounted stripline can potentially be 
exploited in the development of the UHF antenna for AROD. 
[15] A simple antenna to employ is a dipole ‘extending 
from the top of AROD or out the bottom. At a transmission 


frequency of 800 MHz, the wavelength is 


А = E = 14.76 inches (4.1) 


where c = 1.18 х 101% inches/second. 


A half-wave dipole, for example would extend 7.38 inches 
above or below AROD. This design is not considered here. 
Instead, an antenna which conforms to the AROD body is 
developed to satisfy two of the needs listed above: the 
need for a Әбден” fixed antenna and the presence of a 
conducting surface (the body). 

The need for a rugged, fixed antenna is obvious and will 
not be discussed further. The conductivity of the AROD 
body, however, can possibly enhance the conformal antenna's 
radiation pattern characteristics if the body behaves like a 
reflecting ground plane. This potential and basic antenna's 
which can take advantage of the reflecting body are dis- 


cussed in the next section. 


iig 


C. THEORETICAL DEVELOPMENT 

Two basic designs for the UHF antenna are considered 
here. Referring to Figure 4.1, the first is a vertically 
polarized dipole or dipole array and. the е 15 а 
horizontally polarized loop antenna. The dipole array and 
the loop antenna are both very well developed in References 
16 and 17. This work does not endeavor to recreate or 
expand on either books' оо СА development and the 
reader is encouraged to consult these sources if such 
information is desired. Instead, a brief explanation of the 
important characteristics of the dipole and the loop antenna 
is offered as a preface to the practical considerations 
which drive the design process. 

1. Dipole Array 

a. The Half-Wave Dipole 
An omnidirectional far field pattern is the 

goal. For this reason, a dipole which radiates an omni- 
directional pattern is desired. An obvious choice is the 
half-wave dipole. Fed at the center, the current on a half- 


wave dipole will decrease sinusoidally as it nears the end 


points. 

A single half-wave dipole yields the far field 
pattern of Figure 4.2. Тһе horizontal (azimuth) pattern is 
omnidirectional. If the X direction of the dipole cor- 


responds to the X direction of AROD then it is the horizon- 


tal pattern of the dipole which is of greatest interest in 
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Figure 4.2 Far Field Patterns of a Half-Wave Dipole 
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obtaining an acceptable solution for AROD. Therefore, the 
emphasis is placed on obtaining an acceptable horizontal far 
field — from subsequent antenna designs. 
b. Circular Array of Half-Wave Dipoles 

A dipole which is close to the AROD body will 
have radiation characteristics far different from the free 
space or free standing half-wave dipole. The conducting 
body will, at the very least, block the dipole pattern in a 
direction "behind" or toward the body from the dipole. An 


array of dipoles which conforms to and encircles the body 
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will allow a signal to radiate through 360 degrees in 
azimuth (о). Such an array is called a circular or ring 
array [Torp. 352] A circular array of identical, equally 
spaced half-wave dipoles each excited at its center will 
result in a non-directional far field pattern. Figure 4.3 
is the resulting far field gain pattern of a four dipole 
ring array. Although the pattern is plagued by peaks and 
valleys creating a wide dynamic range in directive gain, it 
is a good starting point for the design. 

The theory and analysis for circular arrays is 
well developed and an array yielding an acceptable omni- 
directional pattern can be derived using the procedure given 
in [16:p. 354]. However, no procedure is found for an array 
with an obstructing body in its center. For this reason, a 
model of the AROD body which approximates its electromag- 
netic characteristics is useful. To such a model, various 
arrays can be added and the resulting far field patterns can 
ke analyzed for adequacy. 

2. Loop Antennas 
The loop antenna is an attractive alternative för 
the AROD application. A loop would naturally conform to the 
body circumference and could be constructed easily. Loop 
antennas are thoroughly discussed in Reference 17. However, 
analytical results are readily derived only for "small" or 


"moderate" sized loops. A moderate size loop is defined as 


2Pi 


o de е (2.2) 
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Figure 4.3 Horizontal Far Field Pattern - 4 Dipole Ring 


where 
Pi = 3.1415962, 
b is the radius of the loop, 
A is the wavelength. 
A circumferential loop on AROD has a radius b > 15 inches. 


At a transmission frequency of 800 MHz, the wavelength X = 


14.76 inches and equation 4.2 becomes: 


= 
= = 


e b» 6.39 (f = 800 MHz) (4.3) 
Clearly,  AROD does not satisfy the requirements for a 
moderate size loop. | 

The difficulty with deriving analytical results for 
larger than moderate size loops lies in evaluating the 
Maxwell's Equations. Larger loops such as an antenna around 
the AROD body require a much more extensive mathematical 
development to evaluate the equations than is desired here. 
Therefore, a numerical technique with the aid of a computer 
is best suited for evaluating the radiation pattern of the 
proposed design. 

Figure 4.4 is the far field gain pattern of a single 
loop antenna with a radius b = 15 inches. While this 
pattern is apparently directional and inadequate for the 
AROD application, combined with the AROD body model, a 
better pattern may result. 

3. Ground Plane Effects 
Perfect ground planes are defined as infinite, 


plane, and perfect conductors. Practical approximations to 
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a perfect ground plane are planes much larger than A and 
good conductors on the order of 10/ mhos/meter.[16:p. 239] 
Clearly, a cylindrical surface only slightly larger than 2) 
at its greatest dimension can hardly be considered infinite 
and plane. Nevertheless, the fact that the AROD body does 
conduct suggests that some link to ground plane theory can 
be achieved. 1 
Ground plane theory is the study of reflections. Ап 
ideal dipole oriented parallel to a perfect ground plane has 
an image in the plane which is directed out from the plane. 
The effect is to double the image radiated outward and 
perpendicular to the plane. In the case of AROD, the image 
is ian outward from the body in the horizontal plane. 
Further, Since the body surface is curved, the image will be 
reflected at increasingly obtuse angles as distance from the 
dipole increases. The object of modeling the AROD body is 
to determine if this pattern of reflection can be used to 
enhance the far field patterns of Figure 4.3 and 4.4 and 
obtain an adequate antenna design. 
a. Ground Planes and Phase 
Before proceeding, an important consideration is 
phasing of the signal. The 800 to 900 MHz waveform planned 
for the video signal carries with it a period and phase. 
The phase of this waveform will vary as the distance from 
the antenna increases. Reflection theory stipulates that 


the distance a waveform travels to and from a reflecting 
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surface can be viewed as a straight line distance from the 
original source, in this case the antenna. Thus, a phase 
change is associated with any reflection. This phase change 
can act to complicate or, at 180 degrees, it can Cancel the 
original signal. To minimize any such complications, the 
antenna should be positioned flush to the reflecting surface 
or at “aA distance resulting in zero phase change in the 
reflected signal. Further work, ii considered here, might 
utilize phasing as a means to steer the pattern as desired. 
b. Ground Planes and the Loop Antenna 

One very important result of modelling the body 
as a conductor and approximating it as a ground plane is 
that it will significantly alter the characteristics of the 
loop antenna. In fact, a circumferential loop around the 
cylindrical body might alternatively be considered a line 
source over a ground plane. The "line" is the loop which 
carries a current from its source. Unlike a line source, 
however, this current would be seen on the loop from two 
different directions since the "line" ends are connected 
together. With these important differences in mind, line 
source effects will play a major part in shaping the numer- 


ically generated far field patterns in the next section. 
D. BODY MODELLING AND COMPUTER AIDED DESIGN RESULTS 
The complexity of the body conforming antenna makes use 


of a computer algorithm to model it very practical. 
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Although a simple array of dipoles lends itself to analyti- 
cal solution, a large loop the size of the one proposed does 
not. A computer program written at the Naval Ocean Systems 
Center entitled Numerical Electromagnetics Code (NEC) is 
available on the Naval Postgraduate School's mainframe 
computer [18]. This code is used to obtain the far field 
patterns pictured in this work and is very suitable for this 
study. 


The AROD body, although it can be pictured as a simple, 


conducting cylinder (Figure 4.1), is a very complex electro- 
magnetic structure. One approach to simulate such a 
cylinder is to construct a wire grid outline. Another 
method is to utilize a continuous surface model. Both 


techniques are discussed in [16] and both can be modelled 
with NEC. In fact, because of the symmetry of the AROD body 
and the use of symmetry possible in NEC, the continuous 
surface model is somewhat less cumbersome to create. TRS 
modelling process and the resulting computer code are the 
subject of Appendix D. 
1. Results with Dipole Arrays 

The dipole arrays previously discussed are simulated 
here with the AROD body model. While array effects are 
expected to influence the resulting far field pattern, the 
conducting surface of the AROD will contribute significantly 
as well. All runs utilize a single center feed on a .1 inch 


diameter wire which is divided into 11 equal segments to 
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make a half-wave dipole. The first antenna simulated with 
the body model is the four dipole ring array pictured in 


Figure 4.5. Note that the array radius b = 15.5 inches. 





Figure 4.5 Body Model with 4 Dipole Ring Array 


The half inch beyond the body radius of 15 inches is to 
account for practical considerations of mounting the dipoles 
and to prevent "electrical contact" with the body in the 
computer model. Results of the four dipole simulation are 
shown in Figure 4.6 and should be compared with the free 
standing ring array of Figure 4.3. While the resulting 
pattern is hardly adequate, it does demonstrate the benefi- 
Cial effects of the reflecting AROD body. Directivity has 


been reduced. 


TES 








Figure 4.6 Horizontal Far Field - 4 Diple Ring and Body 


pigure 406 indécates that the body has effectively 
duplicated the pattern of each dipole and doubled the 
radiating horizontal pattern. Since the pattern multiplica- 
tion depends largely on the distance — radiating 
sources, an analysis of the spacing is revealing. The 
circumference of the 30 inch diameter AROD is 94.2 inches. 
The antenna radius b = 15.5 inches yields an array radius of 
97.38 inches. The four dipoles, then, are each 24.3 inches 
apart along the circumference. In wavelengths this distance 
is n.65). This is not the optimal spacing. The spacing 


sought is pictured in Figure 4.7. If the cylindrical 
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Figure 4.7 Pattern Multiplication of Dipoles 
Above a Ground Plane 


surface can be ignored and instead thought of as a plane, 
then the interaction of the dipole patterns can be supposed. 


Figure 4.2 showed the horizontal pattern resulting from a 


single half wave dipole. If the surface reflects, then the 
patterns from two adjacent dipoles will combine. Dipoles a 
half wavelength apart will result in patterns adding exactly 
in phase producing sharp peaks and valleys as seen in Figure 
4.7. Patterns from dipoles a quarter wavelength apart will 
tend to smooth each other. Therefore, spacing close to .25) 
is sought to reduce the dynamic range. Equal spacing is 
needed to maintain a non-directional array as previously 
discussed. 

A circumference of 97.38 inches and wavelength A = 
14.76 inches yields no whole multiple of .25. Table 4.1 


lists the multiples of quarter wavelength spacing 


TABLE 4.1 


QUARTER WAVELENGTH SPACING OF DIPOLES 


Approximate Number of Actual 
Spacing Dipoles Spacing 


Аа 26 3.7 іпсһев 


31/4 9 10.8 inches 


51/4 5. 19.3 inches 
71/4 4 24.3 inches 


circumference= 97.8 inches 
A/4= 3.69 inches (800 Mhz) 





possible on the body. Four dipoles did not yield an 
acceptable pattern. Figure 4.8 is the result of five 


equally spaced half-wave dipoles. The pattern is clearly 
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more smoothed although more complex than the four dipole 
case. The horizontal, vertical and total patterns are shown 
in Figure 4.8, but the total pattern (solid line) is of most 
interest here. In this case the smoothing effects of the 
reflecting body has reduced significantly the sharp drop off 
in directive gain characterized by the four dipole example. 

“Тһе case of a nine dipole array deserves ‘some 
practical discussion about feed points. The problem of 
feeding multiple points in phase from one transmitter, as in 
the case of AROD, is easy to overcome on the computer but 
more difficult with actual components. Finite distances 
exist between transmitter and feed point and phasing is 
Simplified with fewer feed points. Additionally, power 
losses associated with resistive antennas may become too 
great with multiple loads. Therefore, for the nine dipole 
array, the simulation will consider only three equally 
Spaced feed points. 

Figure 4.9 shows the results of this arrangement. 
While not ‘a very omnidirectional pattern, the nine dipole 
case reveals that reducing the number of feed points does 
not increase the directivity and,to the contrary, reduces 
the severity of the peak to valley gain difference of 
dynamic range. With this observation in mind, the five 
dipole antenna which produced a smoother pattern may be 


enhanced with fewer feed points. Figure 4.10(a) is the 


result with only one feed point at » = O on the five dipole 
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array. The far field pattern is greatly improved over the 
five feed case and the dynamic range is greatly reduced. 
Figures 4.10(b) and (c) are vertical cuts over zenith (9) 
for two different azimuth angles. The patterns are non- 
directional and show some promise as a suitable antenna. 
2. Results with Loop Antennas 

Two variants of the loop antenna with AROD body 
model are studied here. As in the case of the dipole 
arrays, reflections play a major role in improving the far 
field pattern over the free space results of Figure 4.4. 

The loop is a single, circumferential wire pictured 
in Figure 4.11. The first variant has one feed point at ф = 
0. The horizontal pattern which results is shown in Figure 
4.12(a) and differs markedly from the free standing loop 
pattern. While the gain peaks of the free standing pattern 
can be detected in Figure 4.12(a) at » = 0, 60, 120, 240 and 
300 degrees, the other peaks which give a desirable non- 
directional far field are not a result of the loop theory 
analysis. However, if the antenna is Жаа аз а ТЕ 
source as previously discussed, this pattern is essential 
an 800 MHz wave propagating along the wire. A peak is seen 
every quarter wavelength with only slight variation at o - 
1809 probably due to opposing currents meeting at this 
point. The currents are not in phase when they meet since 
the distance travelled (half the circumference) is not a 


whole multiple of .5). The 97.38 inch circumference fits 
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Figure 4.10 (b) Vertical Far Field (p=0) 
5 Dipoles with Body · 


132 








Figure 4.11 Single Loop and Body Model 


a 64 line source with a feed at p = 0. The peaks which are 
seen on the far field pattern of Figure 4.12(a) can be 
viewed as the peaks which result from the line source. 

The far field pattern of the single feed loop is 
further described in Figures 4.12(b) and (c). For the case 
when » = 0, the zenith or vertical pattern Бабуна that of 
a long (60) lines source endfire array as described in 
[16:p. 180). Discrepancies between the theoretical endfire 
far field and that pictured here are the result of (1) the 
fact that this antenna is a loop, (2) an actual line source 
lenght of 6.5}, and (3) opposing currents due to the shape 


of the wire as a loop. 
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The second variant is a single loop with two feed 
points fed out of phase. The feed points are at q-0 and 
р=180 degrees. The horizontal pattern is shown in Figure 


4.13 and has more directivity than the single feed case. 


E. ANALYSIS & DESIGN CONCLUSIONS 


1. Analysis 

The goal of the design process was to model the AROD 
body and determine some basic antennas that yield "T— ш 
tional horizontal far field patterns. An iterative approach 
utilizing numerical techniques in the form of NEC was used 
to screen various designs. Only a few of these designs were 
discussed. 

The physical requirement most restrictive to the 
design was that it conform to the cylindrical AROD body. 
This constraint led to the choice of half-wave dipole arrays 
and a loop. The study did not result in a perfectly 
omnidirectional antenna, but two designs which show poten- 
tial for refinement may have been found. | in this 
case is judged by the dynamic range that a radio receiver is 
required to have in order to sense an 800 MHz signal 
regardless of the orientation, in azimuth, of AROD. In 
other words, the horizontal far field pattern of the AROD 


must have a minimum dynamic excursion between peaks. 
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2. Conclusions 

The two designs which present the best potential in 
for refinement are the single loop, single feed and the five 
dipole, single feed array. This conclusion is a result of 
their overall, relatively smooth horizontal pattern. Both 
of these antennas are described in terms of horizontal and 
vertical far field patterns in Figures 4.12 and 4.10. 

The most important assumption which was made, 
relative to the resulting patterns, is that the AROD body is 
a good conductor. This assumption is based on the electri- 
cal characteristics of the carbon fiber composite used to 
construct the body. It should be noted that this conducting 
shield is needed to provide the reflecting surface exploited 
by the antenna designs as well as to shield frequency 

interference between the antenna and the ы” 


internal to AROD. 


— zz, 
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V. CONCLUSIONS 


A. DYNAMIC MODEL 

A computer model of the ducted fan hovering vehicle, 
AROD, was developed by considering the vehicle as both a 
gyroscope and as a flying vehicle. ` Experimental data was 
integrated with the model to simulate the forces and moments 
which act on AROD in flight. Problems with the application 
of some of the data to the computer model were experienced. 
However, the model without the problematic data behaved as 


expected of such a device. 


B. CONTROL 

A linearized approximation of the nonlinear model was 
developed so that optimal control could be used to obtain a 
Steady-state gain matrix with the ultimate goal of controll- 
ing the nonlinear computer model. Optimal control proved a 
viable means of determining a constant gain schedule for the 
multiple-input, multiple-output linear system. The result- 
ing controller also found application to the nonlinear model 


with the problematic data removed. 


C. ANTENNA 
Two basic antenna designs were analyzed for use as an 


omnidirectional antenna for a UHF video downlink on the 
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AROD. A five dipole array and a single loop antenna were 
found to have the best potential of the designs which were 
tested using a computer simulation. The body of the AROD 
was found to significantly alter the directivity charac- 


teristics of the antennas. 


D. FURTHER WORK 


In the areas of dynamic modelling and control of AROD: 


(1) The aerodynamic data and its application should be 
verified. 
(2) A training simulator based on the linear model can 


be implemented utilizing the gain schedule derived 
here on small, readily available computers. 


In the area of a UHF antenna: 


(1) Modelling the entire AROD including the forebody and 
testing the designs for variation from these 
results. | 


(2) Modelling multiple loop antennas with differing 
phasing in an effort to steer the vertical pattern. 

(3) Analyze the antennas for power requirements as 
weight of batteries and power supplies is critical 
to the AROD. 


(4) Test the designs on the AROD prototype for 
suitability. 
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APPENDIX A 


EXPERIMENTAL DATA 
Empirical data received from Sandia National 


Laboratories is listed here. This data was obtained 


experimentally through measurements and wind tunnel testing. 


1. Weights and Moments of Inertia 


Description Name Value(units) 
Weight of Prototype WtPro 85 (lbs) 
Weight of Model WtMod 76.5 (lbs) 
Moments of Inertia: 
about X-axis Ty 1.5246 (ft? lbs) 
about Y-axis - 1.6767 (ft? lbs) 
about Z-axis M. 1.6684 (ft? lbs) 
Products of Inertia:(assumed) I, = І = I = O 
xy XZ yz 
Moment of Inertia of Propeller 
about X-axis 16. „015012 (ft? lbs) 
about Y-axis (assumed) -- 0 
about Z-axis(assumed) -- О 
Propeller Weight (derived) PropWt 4.50017 (lbs) 


2. Hardware Characteristics 


Propeller Efficiency Pref 10.472 (rad/lb-sec) 
Minimum Thrust of Engine Thrmin 35 (lbs) 
Maximum Thrust of Engine Thrmax 14005 (lbs) 
Maximum Throttle Setting TAmax 100 o АРҒЫ 
Maximum Change in Throttle RT max 100 (lbs/sec* ) 
Propeller Speed at Hover HovRat 7200 (rpm) 
Propeller Compensation RHo .0019 (ft-lbx) 

Max Deflection of Vanes MaxDef ¿5236 (tad) 

Max Def Rate of Vanes MaxRatDef .87266 (rad/sec) 
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Description Name Value (units) 

Control Surface Effectiveness _2 
Aileron L eff -21.29 (sec_,) 
Rudder M eff -14.51 (sec_,) 
Elevator N eff -16.68 (sec ) 


Aerodynamic Data 

The data listed in the table below is obtained from wind 
tunnel tests. The manner in which it is applied in the 
simulation (relationship to forces and moments on the AROD 


body) is given in the next section. 


Angle of Attack, AOATOT = V(a* + 6°) 

Name 50 55 60 65 70 75 80 90 

CRdel -- -- -- .0807 .0743 .0712 .881 .085 
Rslope .1514 .1533 .1673 өй -1834 „1764 .2041 .20 
Req 9.293 8.82 8.115 7.848 7.347 7.134 5.011 0.0 
CPdel -- -- -- .3106 .3081 .3305 .3515 0.0 
Pslope .5633 .5123 .6974 .7612 .7603 .8189 .8853 .90 
Ред 43.14 31.00 29.58 26.74 24.55 20.43 15.17 0.0 
CYdel -- -- -- .0024 -.002 .0090 .0120 0.0 
Yslope_.0609 .0663 .0215 .0060 -.004 „0223 .0302 .04 
Yeq „555 .576 | .006  -.09  -1.63 -1.77  -1.74 -1.8 
VanEff -- -- -- 1.5 1.5 1.5 1.9 1.0 
Veq 51.62 44.96 38.56 32.94 27.76) оа вв ив Зони 
CLdel -- -- --  .2428  .2391 .1907 .1734 .1046` 
Lslope .7954 .7183 .6598 .5949 .5901 .4726 .4366 0.0 
Leq 85. 85. 85. 85. 85. 85. 85. 85. 
CDdel -- -- -- .5754 .6052 .5922 .4966. 3890 
Dslope 1.302 1.313 1.371 1.410 1.493 1.467 1.257 1.2 
Deq 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 
CSdel -- -- -- 2.0091 .0258 .0331 .0318 0.0 
XSslope .0801 .0801 .0801 .0801  .0801  .0801 .0801 .0801 
Seq 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 


NOTE: (--) indicates that no data 


under this heading was 
received for this angle. 
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a. Aerodynamic Forces and Moments 


The aerodynamic data listed in paragraph 3, above, 
results in forces moments about the velocity axes, Уо: 
The relationships developed by the engineers who obtained 


the data are reproduced here. Some convenient groupings of 


constants and parameters of the system are given here. 





_ МЕМоа Е 
ЕРЕ терте 


RHOVA = RHo*Pi* (WtRat) *Veq(AOATOT) 

шу = angular speed of propeller (rad/sec) 

UseDel = -ur*RHovA*WtRat 

UseEg = (1 + 2 *ur) *WtRat 

UseSlp - WtRat*(V, , - V(WtRat) *Veg (AOATOT) 
The Forces and Moments which result from ene 


aerodynamic data are given here. 


-FORCES (with respect to Vtot) 


Due to Body Aerodynamics 


F2]= UseDel*CIdel (AOATOT) + UseSlp*Lslope(AOATOT) 
Faq= UseDel*CDdel(AOATOT) + UseSlp*Dslope (AOATOT) 
Fas= UseDel*CSdel(AOATOA) + UseSlp*Sslope 

Due to Thrust Force Aerodynamics 
Fey= Userq*Leq 

Total Aerodynamic Forces 
ee ae) Pe 
Fa= Fad 


Fe = Fas 
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MOMENTS (with respect to V+ot) 
Due to Body Aerodynamics 
Rar= UseDel*CRdel (AOATOT) + UseSlp*Rslope(AOATOT) 
Pap= UseDel*CPdel (AOATOT) + UseSlp*Pslope (AOATOT) 
Уау= UseDel*CYdel (AOATOA) + UseSlp*Yslope (AOATOT) 
Due to Thrust Force Aerodynamics 
Rtr= UseEq*Reg(AOATOT) 
Ptp= UseEq*Peq(AOATOT) 
Yty= UseEq*Yeq(AOATOT) 
Total Aerodynamic Moments 
Rr= Rar + Rtr 
Рр= Рар + Ptp 
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APPENDIX B 


SIMULATIONS SOURCE CODE 


1. Simple Gyroscope Similation 
The following source code is a Dynamic Simulation Language (DSL) 
program ‘to simulate the AROD аз а gyroscope. This simulation is 


discussed in Chapter II. 


TITLE BASIC AROD GYRO SIMULATION 


de 


CONST IRX=69.52,1xXxX=7063.39,1YY=7768.22,122=27729.58,WRxX=784. 
* 
INIT 

PO= 0. 

Q0=0. 

RO= 0. 

MX = 0. 

M2 0. 


* 


DYNAMIC 
IF (TIME .LT. TO) THEN 
MY = 0 
ELSE 
*** STEP INPUT *** 
T1 = TO + DUR 
MY = AMP*(STEP(TO) - STEP(T1)) 
END IF 


DERIV 
НХ = IRX*WRX + IXX*P 
HY = 1YY*Q 
HZ = IZZ*R 
PD (-1.*HZ/IXX)*Q + (HY/IXX)*R + MX/IXX 
ЧО FENZ/ TJ =P - (HX/IYY)*R + MY/IYY 
RD =(-1.*HY/122)*P + CHX/IZZ)*Q + М2/122 


Р = INTGRL(PO,PD) 
а = INTGRL(Q0,QD) 
R = INTGRL(RO,RD) 


* 


PARAM TO=1, DUR=1, AMP=10000 
CONTROL FINTIM=5, DELT=.05 
SAVE .05, P,Q,R,MY 

PRINT .1,MY,P,Q,R 
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GRAPH (GO,DE=TEK618) TIMEC(NI=5,UN=SEC) MYC(UN='LB-IN(S$Q)') 

LABEL (GO) INPUT TORQUE 

GRAPH (G1,DE=TEK618) TIMEC(NI=5,UN=SEC) PCUN='RADIANS/SEC',NI=3,... 
SC=2.,L0=-2.) Q(LI=3,NI1=3,SC=2.,LO=-2.) R(LI=4,NI=3, ... 
5С-2.,10--2.) 

LABEL G2 ROLL, PITCH, YAW RESULTING FROM STEP PITCH INPUT’ 

END 

STOP 


2. Complete Equations of Motion Simulation 

The following source code is the DSL program used to generate the 
results discussed in Chapter II for the complete AROD model. Notation 
follows that of Appendix ADATA wherever possible. Data listed in 


Appendix ADATA is also reproduced here in the form used by DSL. 


TITLE AROD EQUATIONS OF MOTION - FULL MODEL, NO CONTROL - 7 MAY 87 
CONST GRAV= 32.174, Pl= 3.1415962 
ARRAY TRANS(C9), ABTRAN(9) 


* 


* * k k k k k k k k k k k k k k k kk k k k kiki kikiki kikiki ikki kikiki kk kk kk kk kk kk kk kk kk kk k k 


* MNEMONICS * 
* CR/P/YDEL == ROLL/PITCH/YAW COEFFICIENTS * 
* R/P/YSLOPE== ROLL/PITCH/YAW SLOPES OF COEFFICIENT CURVES * 
* THRMIN/MAX== THRUST FORCE LIMITS (LBS) * 
* R/P/YEQ == EQUIVALENT ROLL/PITCH/YAM FORCE FACTOR * 
* VEQ == EQUIVALENT VELOCITY FACTOR * 
* VANEFF == CONTROL SURFACE EFFECTIVENESS ж 
* CL/S/DDEL == LIFT/SIDEFORCE/DRAG COEFFICIENTS * 
* L/S/DSLOPE== LIFT/SIDEFORCE/DRAG SLOPES OF COEFFICIENT CURVES * 
* FRCSLG == CONVERSION FOR IN**2-LB(FORCE) TO FT**2-SLUGS * 
* WTRAT == RATIO OF AROD WEIGHT TO WIND TUNNEL MODEL WEIGHT ж 
ж P/Q/R == ROLL/PITCH/YAW ANGLE RATES FOR EOM (RAD/SEC) «м 
ж UU/VV/WW == BODY-FIXED VELOCITY COMPONENTS (X/Y/2)(FT/SEC) * 
* PHI/THT/PSI== EARTH-FIXED EULER ANGLES (RADIANS) = 
ж ALPHA/BETA/AOATOT== ANGLES OF ATTACK (RADIANS) * 
* DELE/R/A == ELEVATOR/RUDDER/AILERON/THROTTLE POSITIONS BASED ON * 
* IDENTICAL SERVO EQUATIONS (RADIANS) * 
* THROT == CHANGE IN THRUST FORCE BASED ON THROTTLE SERVO (LBS)* 
* ALT == ALTITUDE BASED ON INTEGRATION OF VERTICAL SPEED (FT)* 
* U/V/WERTH == EARTH FIXED VELOCITIES (X/Y/2Z2) (FT/SEC) * 
* GRNSPD/DIST== SPEED AND DISTANCE (NO DIRECTION) OVER GROUND(FT/S)* 
* THRHOV == THRUST FORCE REQUIRED IN HOVER (LBS) * 
* HOVRAT == ANGULAR RATE OF PROPELLER IN HOVER (RAD/SEC) BASED * 
* ON EXPERIMENTAL DATA * 
ж THR == TOTAL THRUST FORCE (LBS) 
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= NEWVEQ == [MPIRICAL DATA FOR VELOCITY CALCULATIONS ui 
* VDELTA == EQUIVALENT VELOCITY FOR FORCE CALCULATIONS (FT/SEC) * 
R RHOVA == WEIGHTING FACTOR FOR AERO COEFFICIENTS BASED ON - 
* VEQ AND AOATOT = 
> DELTIP == CHANGE IN ROTOR ANGULAR VELOCITY (RAD/SEC) " 
> VTIP == ROTOR ANGULAR VELOCITY (DISCRETE) (RAD/SEC) қ 
" USEDEL/SLP/EQ== COMBINATIONS OF VARIABLES WITH WHICH AERODYNAMIC * 
" DATA CAN BE USED FOR MOMENTS AND FORCES (UNITS BASED ON DATA * 
P L/S/DAF == LIFT/SIDE FORCE/DRAG AERO FORCE COMPONENT (LBS) i 
E L Ss D T.F == LIFT/SIDE FORCE/DRAG THRUST FORCE COMPONENT (LBS) 3 
A FAX/Y/2 == AERODYNAMIC FORCES (LBS) Р 
= FTX/Y/Z == THRUST INDUCED FORCES (LBS) E 
= R/P/YAM == ROLL/PITCH/YAW AERODYNAMIC MOMENT COMPONENT (FT-LBS)* 
E R/P/YTM == ROLL/PITCH/YAW THRUST MOMENT COMPONENT (FT-LBS) x 
m LA/MA/NA == ROLL/PITCH/YAW AERODYNAMIC MOMENT (FT-LBS) я 
is LT/MT/NT == ROLL/PITCH/YAW MOMENT DUE TO THRUST/ROTOR (FT-LBS) n 
= GX/Y/2 == BODY-FIXED GRAVITY COMPONENT (FT/SEC**2) Е 
к = 


жак а а F t kt # kt k f k k k t k k k Y k kt t f k k k ft k k f kt f ft ft k f k kk k k f ft kt k kk kik ik k k k k fW k Tr 


aha hr EXPERIMENTALLY OBTAINED DATA ЖЖЖ Х<«ЖжҰ кі «жж k í x í Xx 


* * 


* e f f K KT f K KT fr ft oce ode ode ode ode one oe oe oe ode cde ose oe ode oe oe ode ode ode ose ode ode A AA AA Aaa хх л a a a a a a a a a a a a a a a a a a a a u U * 


* UNITS * 
ж š * 
* MOMENTS OF INERTIA IXX,IYY,IZZ,IRX LB-FT-SEC**2 * 
* FORCES FATX,FATY,FATZ,ETC L8 * 
* MOMENTS LAT,MAT,NAT,ETC LB-FT * 
* ANGLES PSI,THT,PHI I RADIANS * 
* PSIDEG,THTDEG,PHIDEG DEGREES * 
* ANGLE RATES P,Q,R RADIAN / SEC * 
* VELOCITIES UU,VV,WW,UERTH,ETC FT / SEC * 
* Bun. WR | RADIAN / SEC * 
* * 


*de e e de ode ode ode oe ode ode ode ode ode ode ode ode oe ode k f f k f t kt f kk k ff k kt kt kt k kt f k RA AA ARA ARA ARA HA a хх a a fk kk a a a a a a a a k W Y 


* 


CONST IXXM=7063.39,IYYM=7768.22,122M=7729.58,IRXM=69.552, WEIGHT= 76.5 
* MEASURED IN**2-LB АЧ: Ег ЕВ Гое" 2 - LB L8 


* 


CONST PROPWT= 4.50017, PROPEF= 10.472, RHO= .00192 


* LB UNITLESS PROPELLER WEIGHT/EFFICIENCY 
CONST THRMIN= 35., THRMAX= 115. 
* LB LB MIN/MAX THROTTLE POSITION 
CONST TAMAX= 100., RMAXT= 100. 

| * LBS/SEC LB/SEC*SEC MAX CHANGE/RATE OF CHANGE OF THRUST 
CONST MAXDFL= .5236, RMAX= .87266 
* RADIANS RAD/SEC CONTROL VANE MAX DEFLECTION & RATE 


CONST LAEFF= -21.29, MEEFF= -14.51, NREFF= -16.68 

* ROLL, PITCH, YAW VANE EFFECTIVENESS 

CONST H1= 17.77, H2= 157.91, HG1= 219.912, HG2= 24674.126 
*ELEV, RUD, AIL, THROTTLE GAINS/ GYRO GAINS 


* 


CONST WE= 1., K= .S 
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* THROTTLE INTEGRATION LOOP GAINS 
CONST ALPMIN= .8727, ALPMAX= 1.570796, AOAMIN= 50., AOAMAX= 90. 
* ANGLE OF ATTACK LIMITS FOR TABLE DATA REFERENCE 


* 


k********** AERODYNAMIC PARAMETERS FROM WIND TUNNEL TESTS  ****xwxxxx 
* | : 
* NOTE: ALL TABLE DATA IS A FUNCTION OF ANGLE OF ATTACK (TOTAL), 

* AOATOT RANGES FROM 50 TO 90 DEGREES. 


* 


de e de e e de e ode oe oe oe ode cfe e oe ode oe oe e ove ove oe oo ox ROLL DATA Ve e je e e e fe ode eode fe oe e oe e oe oe oe ode fe oe e oe ode oe oe oe * * W * * 


AFGEN  CRDEL 50,.0807, 55,.0807, 60,.0807, 65,.0807, 
70, .0743, 75,.0712, 5804, ловио 90,.0850 
50,.1514, 55,.1533, 60,.1673, 65,.1977, 
70,.1834, 75,.1764, 80,.2041, 90,.2000 
8,2,16, 50,55,60,65,70,75,80,90, 0,10, 
9.293,8.5, 8.82,7-4, 5 tr "MC 
7.347,5.9, 7.134, 5.4, 50011 ¿AIN | 
He de de de W de de de de de * de de kkk kok xk kk kk PITCH DATA kk kk k kkk kkk kkk wkk kk kkk kkk kkk kkk 
AFGEN CPDEL 50,.3106, 55,.3106, 60,.3106, 65,.3106, 
70,.3081, 75,.3305, 80,.3515, 90,.0000 
50,.5633, 55,.5123, 60,.6974, 65,.7612, 
70,.7603, 75,.8189, 80,.8853, 90,.9000 
8,2,16, 50,55,60,65,70,75,80,90, 0,10, 
34.138,11.69, 31.002,10.45, 29.576,9.72, 26.738,9.44,... 
24.549,6.78, 20.427 /4 766, m5 .1687- OIE MOCOS 
Ve 9e de W W W W de W e W W W W W W W W W e e he A he ox YAM DATA ЖЕТ ж e oe ose e W e «ТЛ «ТТ Ж ЖҰ Жж 
AFGEN CYDEL = 50,.0024, 55,.0024, 60,.0024, 65,.0024, ... 
70,-.0015, 75,.0090,° 80,.0120, 90,.0000 
50,.0609, 55,.0663, 60,.0215, 65,.0060, 
70,-.0038, 75,.0223, 80,.0302, 90,.0400 
8,2,16, 50,55,60,65,70,75,80,90, 0,10, 
0.555,0.04, 0.576,-.51, 0.006,-.46, -.09,-.49,... 
-1.634,-1.3, -1.771,-.55, -1.736,-1.1, -1.8,-1.2 
w kw k k k Ye Ye e e e e de se be be e CONTROL SURFACE EFFECTIVENESS k k ye ye e je e e e e je e e e W W e e e Ж 
AFGEN VANEFF s 50,1.5, 55,1.5, 60,1.5, 65,1.5, ... 
70,1.5, 75,15, 805 lao 90,1.0 

Ye de de de de de de de de de de de de de de Ye de EK VELOCITY COMPENSATING FACTOR W ye de ен ои 
FGEN2 VEQ = 8,2,16, 50,55,60,65,70,75,80,90, 0,10, ... 

51.616,44.93, 44.955,39.11, 38.564,33.00, 32.935,27.45,... 

27.757,21.31, 21.829,15.06қ» 154344 , MIES 0. ORO! 
жж «к к ж e W e e W e je de e e e de oe LIFT DATA k k k W e e e de e e e e e se W W e fe W de e W W W e fe e W fe he e xx 
AFGEN CLDEL = 50,.2428, 55,.2428, 60,.2428, 65,.2428, 
70,.2391, 75,.1907, 80,.1734, 90,.1046 
50,.7954, 55,.7183, 60,.6598, 65,.5949, 
70,.5901, 75,.4726, 80,.4366, 90,.0000 
CONST LEQ = 85. 
ж ж de de sje e e W e je e W e e e e e he e * w DRAG DATA Ve e 3e e e je e e e e de je e e je e e je e e e je de se e je e de de oe e W 
АРСЕН CDDEL = 50,.5754, 55,.5754, 60,.5754, 65,.575%, 
70,.6052, 75,.5922, 80,.4966, 90,.3890 
50,1.3017, 55,1.3131, 60,1.3711, 65,1.4103, 
70,1.4933, 75,1.4674, 80,1.2507, 90,1.2000 
0.0 


AFGEN RSLOPE 


FGEN2 REQ 


AFGEN PSLOPE 


FGEN2 PEQ 


AFGEN YSLOPE 


FGEN2 YEQ 


AFGEN LSLOPE 


AFGEN DSLOPE 


CONST DEQ 
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* f k k k k k k k k k k k & k k k f £f ⁄ í í⁄ W X x *f W WW WW W f f. k. NK f k. f. í. í. жж 


SIDE FORCE DATA 


REGEN CSDEL = 50,.0091, 55,.0091, 60,.0091, 65,:0091, 
70,.0258, 75,.0331, 80,.0318, 90,.0000 

CONST SSLOPE = 0.0801 

CONST SEQ = 0.0 


LR K. AAA AAA AAA AAA AAA AAA AAA AAA AAA AAA AAA AAA AAA AAA AAA AAA AAA AAA AAA AAA u 


INITIAL 


* 


***  CONVERSIONS 
***** FORCE TO SLUGS AND IN**2 TO FT**2 
FRCSLG= 1./(GRAV*144.) 
IXX= IXXM*FRCSLG в. 
IYY= IYYM*FRCSLG 
122= IZZM*FRCSLG 
IRX= IRXM*FRCSLG 
***  MOMENTS OF INERTIA 
IYMIZ= IvY - 122 
IZMIX= IZZ - IXX 
IXMIY= IXX - IYY 
* 
MASS= WEIGHT/GRAV 
WIRAT= MEIGHT/85. 
** f W ** W W W í *# f f he e I de de e a n vx INITIAL CONDITIONS Yo oe T T Еи 
PO = 0. 
00 = 0. 
RO = 0. | 
00 = 0. 
vO = 0. 
но = 0. 
PHIO= 0 
THTO= 0 
PSIO= 0. 
DELEDO= 0. 
DELEO = 0. 
DELRDO= 0. 
DELRO = 0. 
DELADO= 0. 
* DELAO = 0. 
DLTRDO= 0. . 
DLTRO = 0. 
THROTO= 0. 
ALTO = 0. 
GRNSPO= 0. 
THRHOV= WEIGHT 
HOVRAT= (THRHOV - PROPWT)*PROPEF 
* * W KRKKKKKKKK DEFINE MACRO TO LIMIT ANGLE Aw N A A A A dn oe on ove Ye Ye 9e o ox ox x o o X 
MACRO ANGOUT= ANLIMCANGIN, ANGMIN, ANGMAX ) 
ANSIGN= SIGN(1.,ANGIN) 
ANGOUT= LIMITCANGMIN, ANGMAX, ANSIGN*ANGIN) 
ANGOUT= ANSIGN*ANGOUT 


ENDMAC 


яяяяхяяяя 


DEF 


INE MACRO FOR USE 
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IN SERVO EQUATION LIMITING 


ye y y y y y y Y + + 


MACRO OUT= SERVOCOUTO,OUTD,LIM) 
OUTs INTGRL(OUTO,OUTD) 
OUT= LIMITC-LIM,LIM,OUT) 
ENDMAC 


de e e e e de ft ft t k f k f k k f f ne ve e ve e e e e de e ne cde de fe e de e on de de e sn ve e f ff f k k k k f £ f £ f £ f fW f f“ f f f f ууу I NN 


DYNAMIC š " 
* 
жеккен кекке COORDINATE TRANSLATION COMPONENTS сити 
***** BODY-FIXED TO EARTH FIXED (EULER ANGLES) 
* ROW 1 
TRANSC1)= COSCPSI)*COSCTHT) 
TRANSC2)= COS(CPSI)*SINCTHT)*SINCPHI) SINCPSI)*COS(CPHI) 
TRANS(3)= COSCPSI)J*SINCTHT)*COSCPHI) + SINCPSI)*SINCPHI) 
* ROW 2 
TRANS(4)= SINCPSI)*COSCTHT) 
TRANSC5S)= SINCPSI)*SINCTHT)*SINCPHI) 
TRANS(6)= SINCPSI)*SINCTHT)*COSCPHI) 
* ROW 3 
TRANS(7)= -SINCTHT) 
TRANS(8)= COSCTHT)*SINCPHI) 
TRANS(9)= COSCTHT)*COSCPHI) 


+ 


COS(PSI)*COSCPHI) 
SINCPHI)*COSCPSI) 


* 
de cede ee de de e VELOCITIES RA dede dede de de de de de de de de de de de de de de de de de de e e e 
УТОТ = SQRTCUU*UU + VV*VV + WW* WW) 
UERTH= UU*TRANS(1) + VV*TRANS(2) + WW*TRANS(3) 
VERTH= UU*TRANSC4) + VV*TRANS(5) + WW*TRANS(6) 
МЕКТН= ЏОХТКАМ5 (7) + VV*TRANS(C8) + WW*"TRANS(CY) 
GRNSPD= SQRTCVERTH*VERTH + WERTH*WERTH) 


* 


wkkwk wkwkwkwkwk kw wk kk kk wrk kkk kr ANGLE OF ATTACK * * * de Ye e e de de e de e ode ode ode oe e oe e oe e e e ve KH 


* IF(VTOT.LE..5) THEN 
ALPHA= 1.570796 
BETA = 0. 
* ELSE 
* ALPHA= ASINCWW/VTOT) 
* BETA = ASINCVV/VTOT) 
* ENDIF 
ALPHA = ANLIMCALPHA,ALPMIN,ALPMAX) 
BETA = ANLIMC(BETA,0.,1.570796) 
АОАТОТ= RADEG*SQRT(CALPHA*ALPHA + BETA*BETA) 


AOATOT= LIMITCAOAMIN,AOAMAX,AOATOT) 
* 
* * k k w. * k t k ANGLE OF ATTACK COORDINATE TRANSFORMATION “%****w*w*wwww 
*** BODY-FIXED TO VEIROCI TY 
* ROW 1 
ABTRAN(1)= -COSCALPHA)*COSCBETA) 
ABTRAN(2)= SINCBETA) 
АВТААН (3) = SINCALPHA)*COSCBETA) 
* ROW 2 
ABTRAN(4)= COSCALPHA)*SINCBETA) 
АВТКАНМ (5 ) = СОЅ (ВЕТА) 
ABTRAN(Ó)= -SINCALPHA)*SINCBETA) 
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* 


ROM 3 


ABTRAN(7)= -SINCALPHA) 
ABTRAN(8)= 0. 
ABTRAN(9)= -COSCALPHA) 


* 


* k f k k k k k k k k k k k k k kf ft жж * w k e oe ode f fW f ух ft tf YW? Y f? e o N ose кой 


FORCE COMPONENTS 


THR = LIMITCTHRMIN, THRMAX, THRHOV+THROT) 
NEWVEQ= FGEN2(VEQ,AOATOT,O.) 
VDELTA= VTOT- SQRTCWTRAT)*NEWVEQ 
RHOVA = RHO*PI*SQRTCWTRAT) *NEWVEQ 
* 
DELTIP= (THR - WEIGHT)*PROPEF Ж 
VTIP = HOVRAT + DELTIP 
* 
USEDEL= -DELTIP*RHOVA*WTRAT 
USESLP= VDELTA*WTRAT 
USEEQ = (1. + (C2./VTIP)*DELTIP.) * WTRAT 
*** FORCES 
LAF = USEDEL*AFGEN(CLDEL,AOATOT) + USESLP*AFGENC(LSLOPE,AOATOT) 
SAF = USEDEL*AFGEN(CSDEL,AOATOT) + USESLP*SSLOPE 
DAF s USEDEL*AFGEN(CDDEL,AOATOT) + USESLP*AFGEN(DSLOPE,AOATOT) 
LTF = USEEQ*LEQ 
* 
РАХ = ABTRANC1)*DAF + ABTRAN(2)*SAF + ABTRAN(3)*LAF 
FAY = ABTRAN(4)*DAF + ABTRAN(5)*SAF + ABTRAN(6)*LAF 
FAZ = ABTRAN(7)*DAF + ABTRAN(8)*SAF + ABTRAN(9)*LAF 
FTX = ABTRANC3)*LTF 
FTY = ABTRAN(6)*LTF 
FT2 - ABTRAN(9)*LTF 
FATX = FAX + FTX 
FATY = FAY + FTY 
FAIZ = РА? + ЕТ? 
***  MOMENTS 
RAM = USEDEL*AFGENCCRDEL,AOATOT) + USESLP*AFGEN(RSLOPE,AOATOT) 
+ USEEQ*FGEN2CREQ,AOATOT,O.) 
PAM s USEDEL*AFGEN(CPDEL,AOATOT) * USESLP*AFGEN(PSLOPE,AOATOT) 
+ USEEQ*FGEN2CPEQ,AOATOT,O.) 
YAM = USEDEL*AFGEN(CYDEL,AOATOT) + USESLP*AFGEN(CYSLOPE,AOATOT) 
+ USEEQ*FGEN2(YEQ,AOATOT,O.) 
* 
LA = ABTRAN(1)*RAM +ABTRANC2)*PAM +ABTRAN(3)*YAM 
MA = ABTRAN(4)*RAM *ABTRAN(5)*PAM *ABTRAN(Ó6)*YAM 
NA = ABTRAN(7)*RAM +ABTRAN(8)*PAM *ABTRAN(9)*YAM 
LT = -DELTIP*DELTIP*FRCSLG*.0729 + IXX*LAEFF*DELA 
MT = MEEFF*IYY*DELE 
NT = NREFF*IZZ*DELR 
LAT = LA + LT 
MAT = MA + MT 
NAT = NA + NT 
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k kk kk kk kk kk k kk k kk kkt kkk kk GRAVITY COMPONENTS *e e RR RARA ARA RAR RARA RR RR 


GX = -GRAV*TRANS(1) 
GY -GRAV*(-TRANS(4)) 
GZ -GRAV*(-TRANS(7)) 


*e e ve k k k e de de de de de de de de e de de de de ve fe de ve de de de de ye de de de de de de de se de de de de de de e de de ye ve de de k de de de ve 9e ve ve He ve A de ve He de he ve ve e 


* oo RR + + + + + + + + + + + + + + + + + + + + + + + B 4 о + + + + + + + + + + + + + + 
SAMPLE 
ze INSERT CONTROLLER HERE 


* 
* 
*++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ ++ +++ +++ 


DERIVATIVE 


— 


3 

хакке какая яая яяя MOMENT DERIVATIVES  # F-X k d d de de ded de dede de de de e e je e e e e e eee 
PD = (IYMIZ*R*Q - IRX*DELTIP + LAT)/IXX 
QD = (IZMIX*P*R - IRX*VTIP*R + МАТ)/ТҮҮ 
RD = (IXMIY*P*Q + IRX*VTIP*Q + NAT)/IZZ 

de e de de e de oe de de de de e je ee de e e eoe FORCE DERIVATIVES  *****dedede dede c de de de de de de de de te ae dee e eee 
UD = (VV*R - WW*Q) + FATX/MASS + GX 
VD = (WW*P - UU*R) + FATY/MASS + GY 
WD = (UU*Q - VV*P) + FATZ/MASS + GZ 

* k k k * k k kk kk kkkkkkAkkAk OANGLE DERIVATIVES “RR RRR RRR e de de dee e e dee 


PHID= P + (Q*SINCPHI) + R*COSCPHI))*TANCTHT) 

THTD= Q*COSCPHI) - R*SINCPHI) 

PSID= (Q*SINCPHI) + R*COSCPHI))/COSCTHT) 
ж k k k k k k k k k k k k e dee kk SERVO DERIVATIVES k k k w k k de k kk k fe de de He de de de de se de de e k k he de be + 
*** ELEVATOR, RUDDER, AILERON, THROTTLE 

DELEDD= -H1*DELED H2*DELE + H2*UE 

DELRDD= -H1*DELRD H2*DELR + H2*UR 

DELADD= -H1*DELAD H2*DELA + H2*UA 

DLTRDD= -H1*DLTRD H2*DLTR + H2*UT 


THROTD= -WE*THROT + K*WE*DLTR 

И | 

*** INTEGRATIONS 

de Je de de e de de de de de de de de de de de k k k k k k kk MOMENT EQUATIONS de *e e de de de de de de de de de de de de de de de de de fe de de de de * 
P- = INTGRL(PO,PD) 
Q = INTGRL(Q0,QD) 
R = INTGRL(RO,RD) 


* k k k k k k К} FORCE EQUATIONS k k k k k k k k de de de de de k k k k k k de de he KH de x 


UU= INTGRL(U0,UD) 
VV= INTGRL(VO,VD) 
WW= INTGRL(WO,WD) 
*de *e de de de de de de de de de k k k k k k k de de k k k k k + EULER ANGLES ж k de de de de de he de de de e de de de de de de de de de de she he de e 
PHI = INTGRLCPHIO,PHID) 
THT = INTGRLCTHTO,THTD) 
PSI = INTGRL(PSIO,PSID) 
PHIDEG= RADEG*PHI 
THTDEG= RADEG*THT 
PSIDEG= RADEG*PSI 
de Y de de de de de de de kk k de de de w k de de he be de he k k k SERVO EQUAT IONS k k k k k k de be de de de be k be k de be k de He ve k e * k 
DELED SERVO(DELEDO,DELEDD,RMAX) 
DELRD SERVO(DELRDO,DELRDD,RMAX) 
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DELAD = SERVO(DELADO,DELADD,RMAX) 
DLTRD = SERVO(DLTRDO,DLTRDD,RMAXT) 
DELE = INTGRL(DELEO,DELED) 
DELR = SERVOCDELRO,DELRD,MAXDFL) 
DELA = INTGRL(DELAO,DELAD) 
DLTR = SERVOCDLTRO,DLTRD,TAMAX) 
THROT = INTGRLCTHROTO,THROTD) 

*** ALTITUDE 
ALT = INTGRLCALTO,UERTH) 


DIST= INTGRL(GRNSPO,GRNSPD) 
* 
CONTRL FINTIM= 5 
PRINT .5, PHIDEG,THTDEG,PSIDEG,P,Q,R,UU,VV,WW,DELA 
SAVE .05, PHIDEG, THTDEG,PSIDEG,P,Q,R,UU,VV,MM 
PARAM UE =0.,UR =0.,UA=-.052,UT= 0.,DELA0=0 
LABEL (G1) TIME RESPONSE TO UA =-.052(NO CONTROL) 
GRAPH (G1,DE=TEK618) TIME(NI=5,UN=SEC) PHIDEG(NI=4,SC=45,L1=1, 
UN='THTDEG AND PSIDEG (DEGREES)',LO=-90) THTDEGCAX=OMIT,NI=4, 
SC=45,L1=3,L0=-90) PSIDEGCAX=OMIT,NI=4,SC=45,L1=4,L0=-90) 
GRAPH (G2,DE=TEK618,0V,P0=0,5) TIME(CNI=5 ,UN=SEC) PCNI=4, 
UN='Q AND R CRAD/SEC)' ,LO=-2,SC=1,L1=1) Q(AX=OMIT, ... 
NI=4,L0=-2,SC=1,LI=3) R(AX=0MIT,NI=4,LO=-2,SC=1,LI=4) 
GRAPH (G2,DE=TEK618,0V,P0=6,0) TIMEC(NI=5,UN=SEC) VV(NI=4, 
UN='UU AND WW (FT/SEC)',LO=-2,SC=5,LI=3) UUCAX=0MIT, 
N1=4,LO=-2,SC=5,LI=1) WWCAX=OMIT,NI=4,LO=-2,SC=5,L1=4) 
END 
STOP 
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APPENDIX C 


CONTROLLER AND LINEAR MODEL SIMULATION 


The DSI/VS source code used to simulate the linearized model is 
listed here. The steady-state controller is implemented under the ` 


"SAMPLE" segment. This same controller is used in the nonlinear model. 


TITLE SIMPLE MODEL WITH OPTIMAL GAINS FROM CONTROLS - 9 AUG 87 

*** FEED BACK GAINS USING THE ERROR STATE CONTROL LAW U= -FC(X-R) ка 

ARRAY RKPT(8),RKQR(8),XPTC(8),XQR(8), UPT(2),UQR(2),EPT(8), EQR(8) 

FGEN2 FBKPT= 2,8,16, 1,2), 2 (See, 5 266 7 NBN 
-1.99484,-.60299,-.71717,.27437,1.74559, .00994, .14049, .00016, 
-.58535,111.16428,-.20533,28.90878, .58170,1.66845, .00962, .13667 

* 

FGEN2 FBKQR= 2,8,16, 1,2,1,2,3,4,5,6,7,8, | 
-1.61288,-1.61027,-.55320,.31985,1.05412,-.20752,.12656,-.00264,... 
1.56829,-1.59987,-.30479,-.56615,.10693,1.23477,.00096,.12962 

* 

CONST RMAXT= 100., TAMAX= 100., MAXDFL= .5236, RMAX= .87266 

CONST A1=.0236133,A2=-3.7381,A3=-6.750686,A4=.4830933 

CONST A5=6.784433,A6=.107891,A7=.893727 

CONST LAE=-21.29,MEE=-14.51,NRE=-16.68 

CONST H1= 17.77, H2= 157.91 

CONST PHDES=.0,ALRDES=.00, THDES=.0,PSDES=.0 

CONST RKPT(3)=0,RKPT(4)=0,RKPT(5)=0,RKPT(6)=0,RKPT(7)=0,RKPT(8)=0 

CONST RKQR(3)=0,RKQR(4)=0,RKQR(5)=0,RKQR(6)=0,RKQR(7)=0,RKQR(8)=0 

INIT 
DELEDO= 
DELRDO= 
DELADO= 
DLTRDO= 
DELEO= 0 
DELRO= 0 
DELAO= 0 
DLTRO= 0 
PO= 0 
Q0- 0 
RO= 0 
THRO= 0 
I= 1 
jJ 7 


о о O O 
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Ecc DEPINE MACRO FOR USE IN SERVO EQUATION LIMITING 


MACRO OUTs SERVO(OUTO,OUTD,LIM) 
OUT= INTGRLC(OUTO,OUTD) 
OUT= LIMITC-LIM,LIM,OUT) 

ENDMAC 

DYNAMIC 

**** INPUT CALCULATIONS 
RKPT(1) = PHDES*STEP(O.) 


ЕКРТ(2) = ALRDES*STEPCO.) 
RKQR(1) = THDES*STEPCO.) 
RKQR(2) = PSDES*STEPCO.) 


XPT(1)= PHIB 
XPT(2)= U 

XPT(3)= P 

XPT(4)= THROT 
XPT(5)= DELA 
XPT(6)= DLTR 
XPT(7)= DELAD 
ХРТ(8) = OLTRO 
ХОВ (1) = THTB 
ХОК (2) = PSIB 
XQR(3)= Q 

XQR(4)= R 

XQR(5)= DELE 
XQR(6)= РЕГЕ 
XQR(7)= DELED 
XQR(8)= DELRD 


* 


SAMPLE 
* 
00 30 I=1,2 
BUPT =0 
BUQR =0 
DO 20 J=1,8 
ЕРТСЈ)= ХРТСЈ) - RKPT(J) 
EQRCJ)= XQRCJ) - RKARCJ) 


FBPT= FGEN2CFBKPT,I,J) 
FBQR= FGEN2(FBKQR,1,J) 
BUPT= BUPT + EPT(J)*FBPT 
BUQR= BUQR + EQR(J)*FBQR 
20 CONTINUE 


UPTCI)= -CBUPT) 
UQR(I)= -(BUQR) 


30 CONTINUE 


UA= ОРТ(1) 
UT= UPT(2) 
UE= UQR(1) 
UR= UQR(2) 
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de f f de de ff ft fff lr t 


DERIV 


* * * 


PHIBD 
UD = 
THTBD 
PSIBD 


P 


АЕТТНКОТ 


Q 
R 


PD= A1*U + A2*THROT + LAE*DELA 
THROTD= 
QD= A3*R + MEE*DELE 
RD= AS*Q + NRE*DELR 
ELEVATOR, 
DELADD= 
DLTRDD= 
DELEDD= 
DELRDD= 


PHIB= 


U 


THTB= 
PSIB= 


THROT 


- ТНВОТ + „оре 


RUDDER, AILERON, THROTTLE 
-HIYDELAD - NNHP"*DETA- H2 NIA 
- НОТ - НЫЕ + Пу 
-H1*DELED - H2e DEINE + НОЕ 
-H1*DELRD - H2"*DEER PEH o SUE 


INTGRL(PHIBO,PHIBD) 
-INTGRL(UO,UD) 

INTGRL(THTBO,THTBD) 

INTGRL(PSIBO,PSIBD) 
P= INTGRL(PO,PD) 


I 


NTGRLCTHRO,THROTD) 


Q= INTGRL(Q0,QD) 
R= INTGRL(RO,RD) 


DELED 
DELRD 
DELAD 
DLTRD 
DELE 
DELR 
DELA 
DLIR 
DES=. 
ZER= 


1 
0. 


SERVO(DELEDO,DELEDD,RMAX) 
SERVO(DELRDO,DELRDD,RMAX) 
SERVO(DELADO,DELADD,RMAX) 
SERVO(DLTRDO,DLTRDD,RMAXT) 
INTGRL(DELEO,DELED) 
SERVO(DELRO,DELRD,MAXDFL) 
INTGRL(DELAO,DELAD) 
SERVO(DLTRO,DLTRD,TAMAX) 


CONTROL FINTIM=5,DELS=.04,DELT=.005- 
PARAM PHIBO=.00,U0=0, THTBO0=0,PSIBO=0 
.04,PHIB,U,THTB,PSIB, ZER,DES 
.2,PHIB,U,THTB,PSIB 

LABEL (G1) REGULATOR GAINS - ERROR STATE 


SAVE 
PRINT 


ZERCAX=0MIT,NI=4,SC=.05,LI=4,LO=-.05) 


GRAPH (G2,DE=TEK618,0V,P0=6,4.5) TIME(NI=S,UN=SEC) THTBC(NI=4, 
UN='AND PSIB R/SEC',LO=-.05,SC=.05,LI=1) PSIBCAX=OMIT, 
МІ-4,10--.05,5С-.05,11-3) ОЕ5(АХ-ОМІТ,МІ-4,10--.05,5С-.05 


ZERCAX=0MIT,NI=4,SC=.05,LI=4,LO=-.05) 


END 
STOP 
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PHI 


(0)= 


ET 
GRAPH (G1,DE=TEK618) TIME(NI=5,UN=SEC) PHIBC(NI=4,SC=.05,LI=1, 
UN='R/SEC AND ALT RATE F/SEC',LO=-.05) UCAX=0MIT,NI=4, 

SC=.5,LI=3,LO=-.5) DESCAX=0MIT,NI=10,SC=.05,LI=4,LO=-.05) . 


, 


LI=4) 


APPENDIX D 


NUMERICAL ELECTROMAGNETICS CODE 


1. Modelling the AROD Body 
The Numerical Electromagnetics Code (NEC) enables the user to 


— 


simulate a solid surface using either the wire grid or the surface patch 
technique. Extensive use of symmetry in the program and the circular 
symmetry of the AROD model allow that the surface patch approach be 


easily implemented. Figure D.1 illustrates the procedure schematically. 





Figure D.1 Schematic for AROD Body Modelling 
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The following computer code is the NEC model of the body which was 


utilized. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Build the 1st column of patches. The column is 9.929 in angular 
width (arc) and 14 inches high. The arc of the column is in the 
XY plane and the height of the column is in the Z direction at a 
radius of 15 inches from the Z axis. 


SP 0, 0, 15., 0., 5., 18.435, 0., 9.92 
SP 0, O, 15., O., 1.5, 5.711, O., 7.44 
SP Ogel0, 155% бе, -1.5 ЮЛ, с, 2 ч 
$Р 0, 0, 15., 0.,-5.,-18.455, 0., 9.92 


Rotate and duplicate the colum. The 1st column is rotated in the 
XY plane and duplicated to form a column that is now 18.849 wide 
and 14 inches high. 


GM. 0, 1, 0:5 50. aa 055 05 


Add patches to top and bottom to form wedge. These patches 
together with the column form a 18.849 section of the cylinder. 


SP 0, 0, 13.95, 1.16, 7., 90., 4.737, 9.26 
SP 0, 0, 13.95, 1.16,-7.,-90., 4.737, 9.26 
SP 05 0, 1496, 990,57., 90.94.737, 7.94 
$P 0, 0, 11:96, .90%@0,-7..290. 4.737, 7 94 
SP 0, 0, 9.470, .785, 7., 90., 4.737, 9.425 
SP 0, 0, 9.470, .785,-7.,-90., 4.737, 9.425 
SP 0, 0, 4.980, .413, 7., 90., 4.737, 9.92 
SP 0, 0, 6.980, :413,-7.,-90., 4.737, 992 


Complete the cylinder. The entire 18.84° wedge is rotated 19 
times to form a 360° cylinder which is missing only a top and 
bottom center piece. These two center pieces are added along with 
each antenna that is to be simulated. cs scales the dimensions to 
meters from inches and ce declares an end to geometry. 


GR 0, 19 


GS 2 
GE 


Specify frequency. The frequency at which the model and antenna 
are to simulated is specified. In this example, f = 800 MHz. 


FR 0,0,0,0,800.,0. 


Store the model. The entire symmetric model is stored ina 
Green's function matrix for recall when the antenna is simulated. 


WG 


162 


2. Antenna Simulation 


A dipole array and a loop antenna are simulated in the NEC code 


which follows. 


a. 


5 Dipole Ring Array 


A five dipole ring array consists of five half-wave dipoles that 


are arranged in a circle. 


(1) 


(2) 


(3) 


(4) 


Read the Green's function. The AROD body model which is stored in 
a Green's function matrix is recalled to start the simulation. 
This matrix also includes the frequency at which the antenna will 
transmit. 


GF 


Specify the dipoles. Wire .1 inch in radius and 7.38 inches long 
is used to model a half-wave dipole at f = 800 MHz. 11 segments 
in each dipole is used so that NEC automatically calculates 
current losses. 5 dipoles are spaced equally around the circle at 
a radius of 15.5 inches. The dipoles longitudinal axis is in the 
Z direction. 


GW 1, 11, 15.5, 0., 3.69, 15.5, 0.,-3.69, .1 

MN ТТ вето, 14.78, 3.659 4.79, 14.74,-3.69, .1 
ИЕ 1], -12. 54, ОМ, 3.59, -12.54, 9.11,-3.69, .1 
gums, wi, - 1655. ОМЗ. 69, -12.54, -9.11,-3.69, .1 
ВЫ а. 29. - ата 3.69, 4.79, -14.74 953,99, .1 


Complete the cylinder. The two center patches which were excluded 
from the symmetric AROD body model are added, the dimensions are 
scaled to meters, and the geometry is ended with no ground plane. 


EE 0009007 0.7 7., 90, 0., 12.57 
OIDO , D. OR, Du, 905-09, 12657 
Ge 2 
GE 0 


Specify the driving source. In this example, a feed point at the 
center of each dipole provides 1.0 volts. 


Loop Antenna 


A single wire loop antenna is added to the AROD body model. The 


loop is flush to the surface of the AROD. 
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(1) 


(2) 


(3) 


(3) 


(4) 


Read the Green's function. The AROD body model which is stored in 
a Green's function matrix is recalled to start the simulation. 
This matrix also includes the frequency at which the antenna will 
transmit. 


GF 


Specify lst wire in the circle. The circle will consist of wire 
arc segments connected together. The lst segment spans 10° and is 
at a radius of 15.5 inches. The wire is in the XZ plane initially 
and is then moved to the XY plane. 

Gardel 15 5 -5..5. , MEN | 

GM 0, 0, 90. „об. „О ee 1 


Create the loop. The first wire is duplicated and rotated 36 
times to form a complete circle. 


GR 1, 36 


Complete the cylinder. The two center patches which were excluded 
from the symmetric AROD body model are added, the dimensions are 
scaled to meters, and the geometry is ended with no ground plane. 


SP 
SP 
GS 
GE 


, 0, 0., U., Г, ое 
, 0, 0.,"09.,-7.:.,;-905 Cae ter 


Specify the driving source. In this example, a feed point at e = 
0% provides 1.0 volts. 


EXA OA 
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